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SUMMARY 


This  study  examines  the  thesis  that  a  part  ot  the  human 
physiological  system  can  be  simulated  by  a  suitably  constructed 
mathematical  model.  The  model  employed  derives  from  a  class  of 
mathematical  programming  methods  that  were  originally  developed 
for  representing  complex  military  and  Industrial  activities  and 
have  recently  been  used  to  represent  Involved  chemical 
equilibria. 

The  motivation  for  this  research  Is  the  long-range  view 
that  a  successful  mathematical  simulation  of  the  human  system 
or  of  human  subsystems  would  provide  an  Important  tool  for 
biological  Investigations.  A  sufficiently  complex  mathematical 
model — that  Is,  a  model  that  embodies  sufficient  chemical  and 
biological  detail  to  represent  a  whole,  functioning  human  system 
or' subsystem — could  be  used  to  explore  biological  hypotheses, 
environmental  stress  reactions,  and  Interplay  of  dependent 
subsystems,  and  could  serve  as  a  pedagogical  tool  or  even  as  an 
aid  to  medical  diagnosis. 

Of  course,  the  foregoing  long-range  view  Is  an  ultimate 
goal.  For  the  moment,  only  the  techniques,  concepts,  and 
characteristics  of  such  a  mathematical  model  are  being  explored. 

This  paper  presents  the  results  of  a  simulation  of  the 
external  respiratory  function.  Respiration,  and  the  consequent 
gas  exchanges  at  the  lung  surfaces^  Involves  mamy  chemical 
reactions  and  a  transformation  of  venous  blood  Into  arterial  blood. 
This  activity  was  chosen  as  a  test  case  to  explore  the  feasibility 
of  constructing  a  mathematical  model  of  a  human  subsystem. 
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Baaing  their  v/ork  on  the  known  physiological  and  chemical 
aspects  of  this  subsystem,  the  authors  have  constructed  a 
mathematical  model  that,  when  solved,  yields  values  describing 
the  major  phenomena  of  the  subsystem.  The  values  of  some  thirty 
different  molecular  species  as  determined  by  the  model  are  in 
excellent  agreement  with  observed  values . 

In  the  present  form  of  the  model,  emphasis  was  placed  on 
the  chemistry  and  thermodynamics  of  the  subsystem  rather  than 
on  the  physical  dynamics,  such  as  flow  and  mixing.  These  time- 
phased  phenomena  are  expressed  only  Implicitly  at  this  stage. 

In  spite  of  these  and  other  approximations,  this  exercise  appears 
to  demonstrate  that  the  mathematical  art  is  now  capable  of 
representing  such  systems  having  large  and  complex  sets  of 
functionally  interrelated  variables,  and  that  there  is  reason 
to  be  optimistic  about  the  possibilities  of  representing  other 
more  complicated  human  subsystems  and  interconnecting  them. 
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Traditlonally,  the  Doox:or  of  Medicine  has  been  trained 
first  in  the  basic  sciences  of  anatomy,  biochemistry, 
physiology,  pathology,  and  pharmacology.  Subsequently,  in  his 
clinical  years  he  has  been  trained  to  Integrate  and  relate  the 
knowledge  obtained  in  these  basic  sciences  of  medicine  as  they 
apply  to  the  individual  patient. 

Since  the  turn  of  the  century,  the  accumulation  of 
knowledge  —  not  only  in  the  clinical  fields  but  also  in  the 
basic  sciences  —  has  become  so  vast  that  there  has  neces¬ 
sarily  developed  a  fragmentation  into  many  specialties  and 
subspecialties. 

The  human  body  may  be  considered  in  one  aspect  as  per¬ 
haps  the  most  complex  chemical  factory  ever  devised.  It  is  a, 
dynamic  chemical  fac±jQEV-±n  which  there  are  no  absolute, 
values.  The  recent  introduction  of  l3otoi)e  techniques  into 
medicine  has  re— emphasized  the  dynamic  state  of  the  human 
body. 

Recognizing  the  dangers  Inherent  in  oversimplification, 
we  might  nevertheless  say  that  in  medicine  we  are  dealing 
functionally  with  highly  complex  systems  of  oxidation— reduction 
reactions  and  with  control  mechanisms  affecting  the  varying 
rates  of  these  reactions.  Of  the  hundreds  of  parameters 
measured  in  medicine,  we  are,  in  fact,  measuring  indices  of 
rates  of  conversions  at  a  given  time.  The  availability  of 
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energy  to  control  these  rates  may  be  a  critical  factor  In  such 
normal  processes  as  growth  and  aging  —  and  In  the  variations 
from  these  normal  processes,  which  we  may  view  as  diseased 
states  or  pathological  conditions. 

New  approaches  must  be  undertaken  to  assist  In  viewing 
the  human  body  as  a  whole.  These  approaches  must  be  applicable 
to  the  Integration  and  evaluation  of  Information  and  of  a  very 
large  number  of  parameters  that  have  been  accumulated  but  not 
necessarily  Interrelated  In  the  hundreds  of  special  fields  In 
which  research  has  been  undertaken  In  the  study  of  the  human 
body. 

It  now  appears  feasible  to  use  advanced  techniques  of 
mathematical  programming  and  computers  as  one  means  of  gaining 
greater  Insight  Into  the  over— all  complexities  of  the  functioning 
of  the  human  body  In  relation  to  Its  anatomical  structure.  Ihe 
present  trial  study  has  demonstrated  the  feasibility  of  this 
approach  as  applied  to  oxygen  utilization. 

Crawford  P.  Sams,  M.D.,  Research  Physician, 
Institute  of  Engineering  Research 
University  of  California 
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role  of  the  United  States  Air  Force  In  making  possible  the 
present  research  on  a.  human  physiological  system.  Starting  In 
19^9^  the  Air  Force  (with  support  from  the  Army  and  Navy)  con~ 
tracted  with  the  National  Academy  of  Sciences  —  National 
Research  Council  to  gather,  and  to  compile  In  handbook  form, 
basic  established  data  In  the  biological  and  medical  sciences. 
This  work  was  administered  under  the  direction  of  the  Aero 
Medical  Laboratory,  Wright  Air  Development  Center  (WADC).  A 
number  of  outstanding  publications  have  resulted  from  this 
Contract.  Two  of  these  were  heavily  relied  on  to  accomplish 
the  mathematical  representation  developed  In  this  study: 
Standard  Values  In  Blood,  Air  Force  Technical  Report  No.  6059# 
July,  1951,  and  Handbook  of  Biological  Data,  WADC  Technical 
Report  56—275,  Dayton,  Ohio,  October,  1956.  In  addition  to 
containing  extensive  and  reliable  data,  these  publications 
are  unique  in  their  rigorous  presentation  of  the  data  on  a 
comparable  and  representative  basis.  This  rigor  Is  essential 
for  the  construction  and  validation  of  mathematical  models. 

Another  significant  Air  Force  research  contribution  that 
was  a  most  Important  aid  In  the  RAND  work  Is  the  publication. 
Handbook  of  Respiratory  Physiology,  compiled  under  the  aegis 


of  the  USAF  School  of  Aviation  Medicine,  Air  University,  in 
195*^.  This  Is  a  scientific  treatise  on  the  various  aspects 
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of  human  respiration,  with  sections  prepared  by  outstanding  ex¬ 
perts  in  the  field  from  within  and  from  without  the  Air  Force. 

The  contributions  of  Dr.  Crawford  P.  Sams  of  the  University 
of  California  to  the  planning  and  execution  of  this  research 
were  made  possible  by  the  support  of  the  Office  of  Civil  and 
Defense  Mobilization  (OCDM)  through  the  Civil  Defense  Research 
Project  at  the  University.  The  cooperation  of  the  University 
of  California  and  OCDM  in  this  arrangement  la  gratefully  ac¬ 
knowledged  . 
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I.  THE  GENERAL  OBJECTIVE 

The  only  way  of  real  advance  in  biology  lies  In 
taking  as  our  starting  point,  not  the  separated 
parts  of  an  organism  and  Its  environment,  but  the 
whole  organism  In  Its  actual  relation  to  environ¬ 
ment,  and  defining  the  parts  and  activities  In 
this  vfhole  In  terns  Implying  their  existing  re¬ 
lationships  to  the  other  parts  and  activities.  — 

J.  B.  S.  Haldane 

Haldane's  formula tlon^^^  of  the  prerequisite  for 
Important  advances  In  biology  was  made  by  that  eminent  British 
biologist  thirty— seven  years  ago.  Until  very  recently.  Its 
Implementation  had  seemed  well  beyond  the  capabilities  of 
human  accomplishments  This  Is  not  to  Imply  that  biologists, 
physiologists,  and  medical  researchers  have  failed  to  develop 
a  knowledge  of  the  micromechanisms  of  the  reactions  occurring 
In  the  separated  parts  of  organisms.  On  the  contrary,  these 
scientists,  working  In  many  countries,  have  amassed  a  re¬ 
markable  amount  of  data  relating  to  reactions  occurring  In 
organisms  and  parts  of  organisms  ranging  from  the  single  cell 
to  complicated  multifunctioning  organs  of  the  human  system. 

The  problem  of  interrelating  and  synthesizing  these  multitudes 
of  dependent  reactions  into  entire  systems  or  even  Into  sub¬ 
systems  representative  of  the  organism  has,  however,  remained 
unsolved . 

The  diagnostician  has  resorted  to  Intuition  In  attempting 
to  synthesize  In  his  mind  the  results  of  tests  of  Isolated 
functions  In  terms  of  the  subsystem  of  particular  Interest. 
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Physiologists  have  had  no  suitable  means  for  relating  a  new 
understanding  of  the  functioning  of  a  part  of  an  organism  to 
larger  parts  or  to  the  system  as  a  whole.  Chemists  developing 
new  drugs  have  been  handicapped  in  obtaining  an  appreciation 
of  the  possible  effects  and  side  effects  of  these  drugs  in 
vivo. 

The  coming  nuclear  and  space  age  presents  an  especially 
difficult  challenge  to  the  biological  sciences  and  to  the 
medical  arts.  The  stresses  that  may  be  placed  on  organisms  by 
the  completely  new  and  stringent  environments  associated  with 
these  technologies  are  largely  unknown.  An  understanding  of 
the  effects  of  these  environments  on  the  whole  organism, 
especially  on  the  hxaman  system,  will  be  difficult  and  slow  if 
not  impossible  to  obtain  within  the  present  state  of  the 
biological  arts  and  sciences.  This  is  partly  because  of  our 
lack  of  knowledge  of  individual  biological  plienomena  and 
partly  because  of  the  present  lack  of  a  technique  allowing 
the  integration  of  a  very  large  number  of  environmental 
effects,  which  may  operate  primarily  through  a  part,  on  the 
whole  system.  We  shall  not  try  to  decide  in  advance  which 
lack  is  more  important.  However,  mathematical  techniques  and 
computing  facilities  are  now  becoming  available  for  Integrating 
Increasingly  complicated  sets  of  components.  It  is  felt  that 
these  techniques  should  be  pushed  to  the  biological  or  mathe¬ 
matical  limit  —  whichever  cones  first  —  in  the  hope  of 
constructing  mathematical  models  for  Increasingly  large 
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portions  of  the  human  system. 

The  Idea  of  building  mathematical  models  of  biological 
systems  Is,  of  course,  not  new.  Many  Investigators  have 
explored  the  possibilities  In  the  past  and  the  field  Is 
attracting  Incz^aslng  attention  at  present .  The  names  of 
lotka,  Rashevsky,  Henderson,  Mlchaells  among  others  come  to 
mind  as  early  Investigators  who  proposed  mathematical  models 
of  biological  systems  or  vdio  pioneered  In  the  quantification 
of  parts  of  such  systems.  Present  Interest  In  the  subject  Is 
attested  by  the  devotion  of  two  complete  Issues  of  Reviews  of 
Modem  Physics  (January  and  April  1959)  to  papers  dealing  with 
the  quantification  of  biological  processes.  Welnrauch  and 
Hetherlngton  have  recently  reviewed  the  Increasing  activity  In 
the  application  of  electronic  computers  to  the  study  of 
biological  problems.  (J.A.M.A.  January  17,  1959  PP*  120-125.) 
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II.  ON  THZ  FEASIBILITY  OF  A  MATHEMATICAL 

MODEL  OF  TK3  IIUIIAN  SYSTEM 

In  con.'^S-lclering  the  forecolns  situation,  the  authors 
believed  that  the  difficult  biological  problem  of  relating  the 
parts  to  the  vfhole,  and  the  v/hole  to  the  environment,  might  be 
well  served  by  employing  some  of  the  newer  mathematical  pro¬ 
gramming  techniques  that  have  been  applied  successfully  In 
representing  other  systems  of  comparable  complexity  and  number 
of  components.  Thus,  linear  and  nonlinear  programming  models 
have  been  successfully  constructed  to  schedule  the  activities 
of  vast  industrial  and  military  enterprises. 

The  £;ener>nl  approach  Is  to  vlew  the  collection  of  organs 
of  the  body  —  the  heart,  lungs,  liver,  skin,  arteries,  veins, 
etc.  —  as  an  Interrelated  and  Interdependent  physical  set  of 
subsystems  each  of  v/hlch  performs  a  v;ell— defined  function  or 
activity  such  as  the  pumping  activity  of  the  heart.  For  a 

-  r" 

particular  organ,  the  level  of  activity  must  be  consistent 
with  certain  Inputs  from  other  organs  and  must  not  exceed  the 
given  capacity  of  the  organ  —  for  example,  the  pumping  ac¬ 
tivity  of  the  heart  Is  limited  by  Its  physical  size  and  the 
availability  of  oxygen  to  heart  muscles,  and  also  by  the 
volume  of  blood  available  from  the  veins .  The  result  of  a 
given  level  of  activity  Is  certain  outputs  to  other  organs. 

If  It  were  valid  to  assume  that  a  change  in  the  activity 
level  of  an  organ  results  In  a  proportional  change  in  Its  In¬ 
puts  and  outputs,  then  the  resulting  mathematical  model  would 


RM-2519 

6 


be  of  the  pure  linear— progranunlng  type.  It  Is  Immediately 
evident,  hovfever,  that  gross  functions  such  as  the  respiratory 
function  or  the  liver  function  do  not  remotely  satisfy  the 
condition  of  proportional  Inputs  and  outputs.  Indeed,  these 
systems  can  function  under  an  almost  endless  variety  of 
complicated  Inputs  producing  a  correspondingly  complicated 
variety  of  outputs.  Nevertheless,  If  we  broaden  our  defini¬ 
tion  of  activity  to  Include  the  forming  of  various  chemical 
species  within  some  gross  function,  then  these  more  elementary 
functions  In  most  respects  satisfy  such  an  assumption  of 
proportionality,  the  well-known  chemical  law  of  combining 
proportions.  In  other  respects,  these  elementary  functions 
are  quite  different.  Thus  the  amounts  of  various  substances 
formed  when  venous  blood  and  air  come  In  contact  through  the 
alveolar  surfaces  of  the  lungs  must  also  satisfy  known  laws 
for  chemical  equilibria.  In  other  cases,  such  as  blood 
circulation,  basic  relationships  drawn  from  the  science  of 
hydrodynamics  play  an  Important  role. 

Thus  a  human— system  model,  while  possessing  many 
characteristics  of  those  used  successfully  to  program  human 
group  activities,  has  other  characteristics  that  result  In 
a  more  complex  structure.  One  aspect  of  the  mathematical 
system  Is  clear.  It  must  Involve  many  variables,  representing 
levels  of  activity  of  certain  organs  and  rates  of  flow  of 
many  substances,  and  these  must  satisfy  a  large  number  of 


RM-2519 

7 


equations  and  Inequalities  expressing  the  limitations  of 
capacity  and  availability  from  other  organs  and  sources. 

It  Is  the  authors'  belief  that  the  successful  development 
of  mathematical  models  capable  of  handling  such  a  multitude  of 
variables,  and  especially  the  more  recent  mating  of  these 
types  of  models  with  chemical  thermodynamics,  provide  en¬ 
couragement  that  the  complexities  of  physiology  may  be  amenable 
to  mathematical  representation  and  solution. 

Developing  a  complete  and  Integrated  human— system 
representation  Is  quite  a  large  order.  In  the  first  place, 
the  size  of  medical  libraries  alone  attests  to  the  large 
volume  of  special  cases  that  such  a  model  would  be  required  In 
some  sense  to  Interpret.  One  might  nevertheless  argue  that, 
with  luck  and  several  years  of  effort  by  a  large  number  of 
qualified  specialists.  It  might  be  possible  to  handle  the  task 
If  It  were  not  for  one  outstanding  fact  —  the  great  gaps  In 
our  knowledge  of  the  human  system.  To  cite  one  obvious 
example,  we  knovf  the  brain  plays  a  major  role  in  the  control 
of  the  physical  processes  of  the  body,  but  the  exact  func¬ 
tional  relationships  that  define  the  control  can  at  best  only 
be  hypothesized. 

Ilovr  then  can  one  presume  to  build  a  mathematical  model  of 
the  vfhole  human  system?  Tlie  only  way  of  doing  this  Is,  as  in 
all  science,  through  the  combined  use  of  fact  and  hypothesis. 
Tlius,  in  areas  where  the  underlying  mechanism  Is  unknown.  It 
Is  generally  true  that  either  X's  or  Y's  hypothesis 


RM-2519 

8 


I 

can  ba  used  to  fill  the  gap.  Incorporation  of  a 
hypothesis  Into  the  model  may  provide  a  useful  way  of  testing 
the  hypothesis  by  comparing  observations  made  under  a  variety 
of  conditions  with  those  predicted  by  the  model  under  similar 
conditions . 

It  Is  the  opinion  of  the  authors  that  a  mathematical 
model  of  the  human  system  can  be  built  largely  on  the  func¬ 
tional  relationships  derived  from  basic  physical  and  chemical 
laws.  Such  a  model  Is  not  limited  by  the  range  of  empirical 
observational  data.  Rather,  the  mathematical  system  i>ermlt8 
the  exploration  of  wide  ranges  of  situations  that  may  repre¬ 
sent  extremes  of  environmental  or  pathological  conditions. 

Even  If  at  first  the  results  of  such  mathematical  explorations 
frequently  do  not  correspond  to  reality,  much  useful  knowl¬ 
edge  can  be  obtained  In  determining  the  reasons  for  the  lack 
of  correspondence. 

The  remainder  of  this  paper  describes  the  results  of 
applying  these  techniques  to  the  human  respiratory  system. 
Other  human  subsystems  should  also  be  amenable  to  representa¬ 
tion  by  the  techniques  described.  Possibly  In  ascending  order 
of  difficulty,  these  might  be  the  following: 

Total  metabolic  system.  Including  digestion, 
nutrition,  and  excretory  functions. 

Synthesis,  storage,  and  conversion  system,  in¬ 
cluding  functions  of  the  liver,  blood— fonnlng 
tissue,  and  muscles. 

Hormonal  regulatory  system. 
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Central  and  pe^lph•^^al  nervous  system. 

Mechano— skeletal  response  system. 

It  Is  quite  possible,  of  course,  that  —  as  the  research 
results  unfold  —  this  ordering  may  be  changed  or  different 
groupings  of  components  may  prove  more  exi>edlent. 

It  Is  not  knovm  whether  the  simultaneous  operation  of 
models  of  several  human  subsystems  In  detail  would  be  beyond 
the  capabilities  of  currently  available  computing  machinery. 
If  each  should  turn  out  to  be  no  bigger  than  the  respiratory 
model,  there  would  be  no  problem.  This,  however,  seems  un¬ 
likely,  Nevertheless,  because  of  some  unique  characteristics 
of  the  chemical  systems  represented,  and  their  thermodynamics, 
and  the  flexibility  of  the  model,  it  appears  to  be  possible 
to  employ  some  tricks  that  may  allow  an  interesting  number  of 
subsystems  to  be  operated  concurrently  with  present  equipment. 
Moreover,  in  many  cases  the  additive  qualities  of  the 
chemistry  can  permit  parts  or  all  of  a  subsystem  to  be  aggre¬ 
gated,  and  yet  allow  the  essential  relationships  with  the 
rest  of  the  subsystem  or  with  other  subsystems  to  be  main¬ 
tained.  This  phenomenon  is  something  like  looking  at  the 
human  system  with  a  microscope.  A  small  field  that  is 
being  examined  can  be  seen  in  great  detail,  while  the  rest  of 
the  system,  still  operating,  is  involved  only  grossly. 
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III.  A  REVIEW  OF  THE  PHYSIOLOGY  AND  CHEMISTRY 
OP  THE  EXTERNAL  RESPIRATORY  SYSTEM 


To  confirm  the  belief  that  complex  physiological  systems 
cein  be  represented  by  mathematical  programming  techniques,  the 
authors  selected  one  Important  subsystem  of  the  human  organism — 
the  external  respiratory  system — ais  an  exercise.*  This  sub¬ 
system  Includes  the  functions  of  breathing,  the  transfer  of 
gases  to  and  from  the  blood,  the  complex  oxygen  and  carbon- 
dioxide  transport  In  the  arterial  and  venous  blood,  the 
reactions  between  blood  plasma  euid  red  cells,  the  reactions 
within  the  red  cells,  the  transfer  of  oxygen  and  carbon  dioxide 
with  the  body  tissues,  and  the  flow  of  blood  through  heart 
action.  An  important  reason  for  selecting  this  subsystem  for 
mathematical  study  Is  that  It  has  been  extensively  Investigated 
and  a  large  amount  of  quantitative  data  are  available  for 
validating  a  model. 

Figure  1  illustrates  in  a  schematic  fashion  the  gross 
physiological  aspects  of  the  external  respiratory  system  In 
man.  The  ultimate  purpose  of  this  system  Is  to  provide  oxygen 
from  the  atmosphere  to  oxidize  carbonaceous  food,  thereby 
generating  heat  and  supplying  energy  for  performing  work,  and 
to  return  to  the  atmosphere  the  COg  (and  some  of  the  water) 
that  Is  formed  by  these  reactions  occurring  In  the  cells  of 


*The  internal  respiratory  system  relates  to  the  metabolic 
processes  occurring  at  the  cellular  level.  These  processes 
are  not  represented  explicitly  In  the  present  exercise. 
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Fig.  I  — Gross  physiological  aspects  of  the  human 
respiratory  system  (  schematic ) 
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body  tissues.  The  exchange  of  gases  between  the  blood  and  the 
external  atmosphere  takes  place  in  the  lungs  by  transfer  across 
the  extremely  thin  walls  of  the  capillaries  and  the  pulmonary 
alveoli . 

Tlie  air  passages  to  and  from  these  exchange  membranes 
Include  the  nasal  cavities,  the  pharynx,  larynx,  trachea, 
bronchi,  and  bronchioles.  The  lung  tissue  proper  contains 
an  Immense  number  of  Irregularly  shaped  air  spaces,  called 
alveolar  sacs.  The  terminal  twigs  of  the  bronchial  tree  open 
Into  these  alveolar  sacs;  see  Fig.  2.  The  number  of  pulmonary 

g 

alveoli  In  the  lungs  of  man  has  been  estimated  at  7*5  x  10 
The  tremendous  surface  area  presented  by  the  alveoli  for  gas 
exchange  results  In  close— to— equilibrium  conditions  between 
arterial  blood  and  sac— gas  composition  except  under  transient 
stress  conditions. 

The  lungs  play  a  passive  role  In  Inspiration  and  expira¬ 
tion.  Their  change  In  volume  during  this  cycle  is  brought 
about  through  change  In  the  capacity  of  the  thoracic  cavity. 

The  change  in  capacity  of  this  cavity  Is  caused  by  the  descent 
and  ascent  of  the  diaphragm,  and  by  the  movement  of  the  rib 
cage . 

Only  about  one— half  liter,  or  one— tenth  of  the  total  air 
capacity  of  the  lungs.  Is  Inspired  and  expired  during  ordinary 
quiet  breathing.  This  Is  called  the  tidal  air.  Even  after 
the  most  forceful  expiration,  from  1000  to  1500  cc  of  air 
remains  In  the  lungs.  The  space  enclosed  by  the  bronchi. 
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bronchioles,  trachea,  larynx,  pharynx,  and  nasal  cavities  is 
called  the  dead-air  space.  It  has  a  capacity  of  about  I50  cc. 
The  air  in  the  air  sacs  is  termed  the  alveolar  air.  Various 
techniques  have  been  developed  for  measuring  or  estimating  the 
composition  of  this  sac  alrj  which  is  in  equilibrium  with 
arterial  blood  and  which  is  different  from  the  external 
atmosphere . 

The  foregoing  characteristics  of  the  lungs,  especially 
the  small  tidal  vol\ime  in  respect  to  total  capacity  and  in 
respect  to  the  large  residual  vol\ime,  permit  the  lungs  to 
accomplish  the  remarkable  function  of  damping  out  the  large 
fluctuations  of  sac-air  composition  that  otherwise,  during 
the  Intermittent  inspiration— expiration  cycle,  would  subject 
the  flowing  blood  to  large  surges  of  oxygen  (and  alternately  of 
carbon  dioxide).  To  a  first  approximation,  this  lung  charac¬ 
teristic  permits  the  breathing  mechanisms  J:o  be  represented 
as  a  large  volume  of  air  of  sac  composition  in  equilibrium 
with  the  blood. 

Actually,  the  breathing  mechanism  and  changes  of  air 
composition  from  atmosphere  to  lung  sac  are  quite  complex  in 
detail.  The  blood  gives  off  carbon  dioxide  (and  usually 
water)  to  the  lung  sacs  and  picks  up  oxygen.  Thus,  the  air 
leaving  the  sacs  Is  of  a  different  composition  from  the  in¬ 
coming  air.  On  top  of  this  change,  there  are  other  factors 
such  as  diffusion,  mixing,  the  dead— air  space,  the  alternate 
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Inspiration-expiration  cycle,  and  the  blood— flow  rate  that 
all  affect  the  composition  of  sac  air  In  relation  to  atmos¬ 
pheric  air. 

One's  breathing  rate  can  be  controlled  voluntarily  within 
limits,  but  It  Is  essentially  an  Involuntary  act  controlled 
automatically  by  the  rhythmical  discharge  of  Impulses  from  a 
group  of  nerve  cells  In  the  medulla  oblongata  of  the  brain. 
This  respiration  center  Is  Influenced  largely  by  changes  of 
carbon-dioxide  concentration  In  the  blood.  Changes  In  oxygen 
concentration  can  also  Influence  the  breathing  rate,  but  this 
Is  thought  to  be  an  Indirect  Influence  through  chemoreceptors 
lying  close  to  the  heart.  These  chemoreceptors  also  api>ear 
to  be  Influenced  within  limits  by  hydrogenr-lon  and  carbon- 
dioxide  concentration.  In  addition  to  affecting  heart— action 
and  blood— flow  rates  and  Indirectly  Influencing  breathing. 

The  gaseous  exchanges  between  the  blood  and  the  tissues 
are  the  reverse  of  those  taking  place  In  the  lungs.  The 
oxygen  concentration  Is  lower  and  the  carbon— dioxide  concen¬ 
tration  Is  higher  In  the  tissues  than  in  arterial  blood. 
Therefore,  In  passing  through  the  systemic  capillaries,  the 
blood  gives  up  oxygen  and  absorbs  carbon  dioxide.  The  compo¬ 
sition  of  the  blood  returning  to  the  lungs  varies,  depending 
on  the  part  of  the  body  It  has  served.  There  Is  also  some 
admixture  of  arterial  blood  with  the  venous  as  It  enters  the 
lungs.  As  a  result,  the  blood  entering  the  lungs  —  which  Is 
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called  mixed— venous  blood  —  is  higher  in  oxygen  and  lower  In 
carbon  dioxide  than  the  most  degraded  venous  blood,  such  as 
that  coming  from  an  extremity  and  sampled  there. 

The  chemical  phenomena  of  the  external  respiratory  system 
in  which  we  are  especially  interested  may  be  represented  sche¬ 
matically  as  in  Fig.  Here  the  system  is  considered  as 
consisting  of  three  compartments,  corresponding  to  the  lung 
sacs,  the  blood  plasma,  and  the  red  blood  cells.  Following 
first  the  path  of  oxygen  from  the  lung  sacs,  we  note  that  in— 
asmvich  as  the  concentration  of  oxygen  dissolved  in  mixed- 
venous  plasma  is  lower  than  that  dictated  by  the  partial 
pressure  of  oxygen  in  the  lung  sacs,  oxygen  transfers  from  the 
sacs  to  plasma,  increasing  its  concentration  there  and  In  red- 
cell  solution.  This  increased  oxygen  concentration  causes 
the  reduced  hemoglobin  in  the  red  cells  to  combine  with  oxygen 
until  a  new  equilibrium  is  reached  at  the  higher  oxygen  con¬ 
centration'.  Hemoglobin  accounts  for  the  large  percentage  of 
oxygen  carried  by  the  blood;  however,  as  is  also  true  of 
carbon  dioxide,  its  solution  in  plasma  and  red  cells  serves 
as  the  Important  pathway  in  and  out  of  the  blood. 

Conversely  to  the  case  of  oxygen,  the  concentration  of 
carbon  dioxide  is  higher  in  mixed- venous  plasma  than  that  dic¬ 
tated  by  the  partial  pressure  of  carbon  dioxide  in  the  lung 
sacs,  so  carbon  dioxide  transfers  from  the  plasma  to  the  sacs 
until  a  new  equilibrivun  is  reached  at  a  lower  carbon-dioxide 
concentration.  But  carbon  dioxide  is  stored  by  the  blood  In 
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both  plasma  and  red  cells  largely  as  bicarbonate  Ion  (and  In 
association  with  reduced  hemoglobin  In  the  carbamlno  form). 

As  the  carbon-dioxide  concentration  Is  reduced,  the  bicarbonate— 
Ion  concentration  Is  also  reduced.  This  reaction  In  the  red 
cells  Is  catalyzed  by  the  enzyme  carbonic  anhydrase.  As  the 
bicarbonate  concentration  In  the  red  cells  decreases,  an  Im¬ 
portant  reaction  called  the  "chloride  shift"  takes  place 
across  the  cell  membrane;  In  this  reaction,  bicarbonate  Ion 
from  the  plasma  transfers  to  the  red  cells,  and  chloride  Ion 
from  the  red  cells  transfers  to  the  plasma.  The  reverse  of 
this  shift  occurs  at  the  tissues  when  carbon— dioxide  concen¬ 
tration  Increases  In  the  blood.  Here,  bicarbonate  Ion  formed 
In  the  red  cells  shifts  to  the  plasma  and  chloride  Ion  shifts 
to  the  red  cells.  This  shift  results  In  an  Increase  In  the 
carbon— dioxide  carrying  power  of  the  blood.  The  remarkable 
Ion— selective  characteristics  of  the  red— cell  membrane  and  the 
resulting  Gibbs— Donnan  equlllbrliom  make  this  Important  shift 
possible  - 

As  mentioned,  carbon  dioxide  Is  also  carried  by  hemo¬ 
globin.  This  hemoglobin— linked  carbon  dioxide,  while  repre¬ 
senting  only  a  small  part  of  the  total  carried  by  the  blood, 
provides  a  significant  part  of  the  carbon  dioxide  transferred 
from  tissues  to  lung  sacs. 

The  physiology  and  the  chemistry  of  the  external  respira¬ 
tory  system  are,  of  course,  much  more  complex  In  detail  than 
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could  possibly  be  indicated  in  this  brief  description.* 

This  introduction,  hov/ever,  should  provide  the  background 
necessary  for  an  understanding  of  the  mathematical  represen¬ 
tation  of  the  more  important  features  of  this  system. 

An  excellent  detailed  description  of  the  external  res¬ 
piratory  system  is  contained  in  Ref.  4. 
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IV.  A  MATHEMATICAL  MODEL  OF  THE  CHEMISTRY 

OF  THE  EXTERNAL  RESFIR/ITORY  SYSTEJ^i 

DESCRIPTION  OF  THE  MODEL 

From  the  physiological  data  available  and  the  working  of 
the  respiratory  system  as  outlined  above  In  Sec.  Ill,  a 
mathematical  model  was  set  up  to  represent  the  more  Important 
of  the  known  Interrelated  physiological  functions  and  chemical 
reactions  Involved  In  the  human  respiratory  system.  Actually, 
two  models  were  formulated.  Model  I  uses  a  simplified  concept 
of  the  hemoglobin  molecule,  but  otherr/lse  Is  sufficiently 
detailed  for  the  present  purpose.  Model  II,  described  In 
Appendix  D,  Is  based  on  Linus  Pauling's  theory  of  the  complex 
hemoglobin  molecular  structure. 

Figure  4,  v/hlch  Is  a  schematic  Illustration  of  the 
external  respiratory  system,  demonstrates  the  relationship  of 
the  Inputs  and  outputs  of  the  mathematical  model  to  the  actual 
system.  The  a's  and  v's  In  the  figure  refer  to  the  Input 
"elements"  coming  to  the  lungs  In  unit  time  from  mixed— venous 
blood  and  from  the  air.  The  x's  refer  to  the  resulting 
numerous  molecular  species  In  the  arterial  blood  and  In  the 
air  of  the  lung  sacs  as  determined  by  the  solution  of  the 
mathematical  model.  At  the  present  stage,  the  cell  outputs 
are  Introduced  Into  the  model  In  terms  of  the  composition  of 
the  mixed— venous  blood. 
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Cell  input 


F'9  4  The  relolionship  of  the  inputs  and  outputs  of 
Model  1  to  the  respiratory  systen. 
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Model  I  vms  constructed  to  provide  an  accounting  of  the 
mass  of  all  of  the  "elements"  Involved.  The  concentration  of 
each  of  the  molecular  species  (dissolved  compounds,  gases. 
Ions)  vras  established  In  terms  of  Its  mole  fraction.  Since 
the  equilibrium  constants  for  the  many  reactions  were  avail¬ 
able,  It  vras  possible  to  Interrelate,  within  the  model,  the 
mutual  effects  of  changes  In  concentration  of  any  of  the 
species.  Further,  from  the  equilibrium  constants  and  mole 
fractions.  It  was  possible,  on  the  assumption  of  Ideal  solu¬ 
tion  behavior  In  the  dilute  solutions  Involved,  to  establish 
the  thermodynamics  of  the  system  within  the  model. 

Some  previous  studies  at  RAND^^^  provided  a  technique  for 
the  Introduction  of  these  energetics  of  the  system  Into  the 
computational  model  In  an  efficient  manner.  This  method  makes 
use  of  the  Gibbs'  chemical-potential  concept  with  the 
minimization  of  total  potential  energy  of  the  system  defining 
equilibrium,  rather  than  use  of  the  equillbinim  constants  and 
their  more  conventional  relation  to  free  energy.  The  two 
concepts  are  mathematically  equivalent.  .In  the  first  part  of 
Appendix  A,  the  interested  reader  v;ill  find  a  review  of  the 
basic  chemical  and  thermodynamic  concepts  employed  In  the  '  - 
present  development.  ''V  ’  * 

J  • 

The  model  constructed  to  represent  this  system  is  shown  In 
Table  1.  Briefly,  the  input  data  of  the  model  consist  of 
given  quantities  of  the  "elements"  of  mixed— venous  blood  and 
air,  shown  on  the  left  in  the  table.  The  output  data  are  the 
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evaluations  of  the  quantities  Xj  representing  the  amounts  of 
various  molecular  species  found  in  three  compartments  rep¬ 
resenting  lung— sac  air,  blood  plasma,  and  I'ed  cells.  Tlie 
numbers  appearing  in  the  various  columns  for  the  individual 
molecular  species  characterize  the  species  and  are  therefore 
Independent  of  their  quantity.  In  greater  detail,  the  model 
consists  of  the  following  portions i 

1.  An  "Input  Section"  consisting  of  given  amounts  v^^ 

of  essential  "elements"  "i"  of  the  blood  corresponding  to 
typical  mixed— venous  blood  values  (measured  in  moles),  plus 
amounts  of  fresh  air  "elements"  "l,"  giving  total  Inputs 

bj^  “  '^l*  simple  chemical  reactions,  the  "elements" 

could  be  defined  as  the  atomic  types  composing  the  molecular 
species.  It  v/as  more  convenient  for  the  present  application 
to  compose  the  various  possible  molecular  species  appearing 
in  blood  from  a  set  of  "elements"  consisting  of  either  atoms 
or  "groups  of  atoms  in  fixed  proportions,"  such  as  or 
0H“.  ) 

2.  Compartment  I,  labeled  "Air  Out,"  with  columns 
corresponding  to  different  species  of  lung— sac  gases  in 
equilibrium  with  the  arterial  plasma  and  red  cells. 

5-  Compartment  II,  labeled  "Arterial  Plasma,"  with 
columns  corresponding  to  chemical  species  appropriate  to  the 
arterial  plasma.  When  the  model  is  "solved,"  the  corresponding 
values  are  the  equilibrium  amounts  (in  moles)  of  the  species 
In  lung  sacs,  plasma,  and  red  cells. 
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4.  Compartment  III,  labeled  "Red  Cells,"  with  columns 
corresponding  to  chemical  species  contained  In  red  cells.  As 
noted  earlier.  Model  I  uses  a  simplified  structure.  A  more 
complete  theory  Is  contained  In  Model  II;  see  Appendix  D. 

The  mathematical  problem  for  Model  I  la  the  following: 

For  given  values  b^  (l  -  1,  2,  ...,  12),  and  for  the 
detached  coefficient  values  a^j  given  In  Table  1,  de¬ 
termine  equilibrium  amounts  Xj  (j  -  1,  2,  ...,  30),  and 
minimum  free  energy  z  defined  by 

*  III  ' 

where  S  denotes  the  sum  over  all  species  J  In  compart— 
k 

ment  k  (k  =  I,  II,  or  III),  Cj  Is  proportional  to  the 
free  energy  per  mole  of  the  Jth  species,  x,.  =  Z  x.  Is 

k  J 

the  total  moles  In  compartment  k,  and  In  means  natural 
logarithm,  subject  to  the  constraints 

Xj  >  0  ( J  =  1,  2,  . . . ,  30), 

30 

^1  ""  ^1J*J  2,  ...,  12). 

Thus,  the  computational  problem  as  presented  In  the  model 
Is  to  determine  the  amounts  of  each  of  the  chemical  species 
present  In  the  lung  sacs.  In  arterial  plasma,  and  In  red  cells 
at  equilibrium. 
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In  the  paat,  the  determination  of  equilibrium  compo¬ 
sitions  of  such  large  chemical  systems  Involved  computational 
problems  that  have  been  tedious  and  time-consuming}  they  have 
often  been  "solved"  by  assuming  certain  species  to  be  dominant 
in  the  final  solution. In  recent  years,  however, 
several  conceptual  and  computational  developments  plus  the 
availability  of  lax^e  electronic  computing  machines  have  made 
it  feasible  to  attempt  the  solution  of  such  large-scale 

problems  repeatedly,  and  under  any  desired  variation  of  the 

f 

basic  Inputs  into  the  system.'  ' 

THE  FREE-ENERQY  FUNCTION 

The  method  here  used  for  solving  the  system  is  described 
in  detail  in  Ref.  5*  It  consists  of  the  direct  minimization  of 
the  nonlinear  free-energy  function  subject  to  certain  linear 
constraints  Instead  of  the  usual  nonlinear  equlllblrum 
equations  that  express  the  mass-action  laws.  Mathematically, 
the  free  energy  F(X)  of  a  mixture  of  n  chemical  species 
relative  to  some  standard  state  containing  Xj  moles  of  the  J— th 
species  can  be  expressed  as 

F(X)  -  Z  x.F  ,  (1) 

J-1  J 

where  X  «*  (x^^,  Xg,  ...,  x^^)  Is  the  set  of  mole  numbers, 

Fj  represents  the  free  energy  per  mole  of  the  J— th  species 
at  the  given  temperature  and  concentration,  Xj/  x;  these  are 
defined  by 
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— 

P.  -  RT  (c  +  In  ^  )  , 

^  X 


_  u 

X  «  Z  X. , 
1 


•J  "  RT 


+  In  P 


(2) 

(5) 

(4) 


P  »  pressure  In  atmospheres, 

T  -  absolute  temperature  Kelvin, 

R  ■  the  gas  constant, 

P?  •»  free  energy  per  mole  of  the 
^  J— th  species  in  its  standard 

state. 

(For  species  in  liquid  state,  the  pressure  terra  In  P  defining 
Cj  is  omitted  for  P's  near  atmospheric;  and  for  pure  species 
in  the  solid  state,  both  this  term  and  the  term  ln(xj/  x)  are 
dropped.)  - 

Where  there  are  several  compartments  instead  of  one,  as 
in  the  model,  it  is  necessary  to  modify  the  free— energy 
equation  so  that  concentration  is  defined  separately  for  each 
compartment.  The  necessary  change  in  F(x)  is  shotm  in  the 
foregoing  formulation  of  the  mathematical  problem  for  Model  I. 

Real  gases  and  solutions  do  not  conform,  over  the  entire 
range  of  composition  and  of  pressure,  to  the  ideal  situation 
described  by  these  equations.  To  retain  the  advantages  of 
the  simple  form  of  the  equations,  it  has  been  expedient  to 
define  a  thermodynamic  function  called  activity,  a^,  which  Is 


proportional  to  concentration  in  mole  fraction: 
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The  constant  of  proportionality,  y^,  called  the  activity 
coefficient,  approaches  1  as  x^/  x  tends  to  0.  In  this  dis¬ 
cussion,  mole  fraction  x^/  x  is  used  instead  of  activity  ay 
with  the  understanding  that  an  approximation  Is  thereby  Intro¬ 
duced.  The  use  of  mole  fraction  Instead  of  activity  In  the 
present  model  with  the  very  dilute  solutions  and  low  pressure 
Involved  does  not  appear  to  Introduce  serious  error. 

The  free-energy  function  Is  thus  defined,  once  the 
numerical  values  of  the  Cj  are  determined  for  the  various 
species.  In  the  case  of  complex  mixtures  of  gases,  extensive 
tables  are  available  (see,  e.g..  Ref.  8)  containing  the  free 
energy  values  Fj/(RT)  for  a  large  number  of  chemical  species 
and  over  a  large  temperature  range.  Unfortunately,  these  data 
do  not  Include  the  values  for  many  of  the  chemical  species 
present  In  physiological  systems.  There  are,  however,  equi¬ 
librium  constants  available  In  the  literature  for  most  of  the 
Important  reactions  occurring  In  the  external  respiratory 
system. 

Appendix  A  describes  In  detail  the  conversion  of  the 
equilibrium  constants  to  the  c^  values  shovm  In  Table  1. 
VJhenever  possible,  data  xrere  dravm  from  basic  physical  chemical 
tables  rather  than  physiological  sources  to  ensure  their 
applicability  under  various  conditions  deviating  from  those 

of  normal  Individuals. 
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THE  MASS-BALANCE  AND  CHARGE  CONSTRAINTS 

The  determination  of  the  equilibrium  composition  Is. 
equivalent  to  finding  the  nonnegative  set  of  values  Xj  that 
minimizes  F(x)  as  defined  by  (l)  above  subject  to  several 
sets  of  linear  constraints.  Hie  fii’st  of  these  are  the  mass- 
balance  constraints, 

n 

^  ^IJ  ^J  “  ^1 

J 

where  there  are  m  different  types  of  "elements";  a^j  indicates 
the  number  of  units  of  element  1  in  a  molecular  species  J, 
and  bj^  Is  the  total  number  of  element  v/elghts  of  element  1 
present  In  the  original  inputs.  The  values  of  the  aj^j  are  the 
entries  In  the  column  below  each  species  in  Table  1. 

Chemically  this  is  equivalent  to  stating  that  the  moss  of  the 
numerous  chemical  species  must  be  in  balance  with  those  of 
the  Inputs.  For  example,  the  first  mass— balance  equation  can 
be  obtained  from  Table  1  by  multiplying  the  unknowns  Xj  by 
the  corresponding  entries  in  the  "Og"  rov^  to  yield 

aj  +  Vj  -  .  Xj  +  Xj  +  Xj.^  +  Xjg  . 

In  addition  to  the  mass— balance  consti’aints  of  the  model, 
there  is  one  expressed  by  the  "charge"  rov;,  which  Indicates 
the  conservation  of  electrical  charge.  It  is  necessary  to 
introduce  this  Into  tlie  model  to  account  for  the  remarkable 
Ion— selective  characteristics  of  the  membrane  separating  tlie 
plasma  and  red  cell  Is.  The  formulation  of  this  set  of 
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constraints  Implies  the  Gibbs— Donnan  equilibrium  that  occurs 
with  such  Ion— selective  membranes;  It  also  provides  for  the 
so-called  "chloride  shift"  that  occurs  through  this  membrane. 
This  constraint  may  be  formulated  In  the  present  context  as 

(^) 

II  J  J 

where  Zj,  the  constant  appearing  In  the  charge  row.  Is  the 
net  number  of  positive  or  negative  charges  present  per  unit 
of  the  charged  species.  The  summation  extends  over  all 
species  In  Compartment  II. 

In  the  operations  of  the  model  so  far,  the  assumption  has 
been  made  that  there  Is  no  net  charge  difference  In  the  Input 
venous  blood,  and  the  charged  Input  "elements"  were  carefully 
balanced  In  accordance  with  this  assumption.  If,  as  Is  shown 
In  Table  1,  one  of  the  compartments  with  Ions  contains  the 
constraint  that  all  charges  In  the  compartment  sum  to  zero, 
and  If  the  Input  charges  sum  to  zero.  It  Is  easy  to  show  that 
the  mass— balance  equations  Imply  the  electrical  neutrality  of 
the  other  two  compartments.  Since  the  first  compartment  con¬ 
tains  only  electrically  neutral  species,  (6)  Implies  that  all 
compartments  are  electrically  neutral.  Besides  providing  for 
the  representation  of  the  Gibbs— Donnan  equilibrium  across  the 
cell  membrane  and  the  resulting  chloride  shift,  the  charge 

equation  could  be  used  as  a  means  for  forcing  other  compo— 

K 

i* 

sltlonal  changes  betv/een  compartments  v/hen  the  charge 
difference  Is  not  zero. 
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Some  JuOfnnent  entered  into  the  choice  of  the  appropriate 
input  "elements"  and  possible  product  species  present  in 
blood.  As  for  the  latter,  those  wei’o  selected  to  yield  a  fair 
representation  of  the  impoi’tant  chemical  phenomena  that  have 
been  observed  in  the  external  respiratory  system.  Experience 
has  shown  that  the  model  illustrated  in  Table  1  is  about  the 
smallest  matrix  in  terms  of  I’ov/s,  columns,  and  compartments 
required  to  represent  those  moi’e  important  phenomena. 

There  are  also  many  possible  choices  of  "elements"  to  be 
used  to  fomi  the  molecular  species.  The  one  selected  here 
happens  to  be  minimum  in  number  (notice  that  there  are  the 
same  number  of  vectorally  independent  columns  as  there  are 
restricting  equations).  Tlie  selection  was  also  made  so  that 
there  v/as  a  physical  interpretation,  in  tenns  of  nonover— 
lapping  groups  of  atoms,  that  ponnitted  obvious  simple  whole- 
number  combinations  to  form  the  various  species. 
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V.  SOLUTIONS  ON  THE  DIGITAL  COMPUTER 

AND 

SIMULATION  ON  THE  ANALOG  COMPUTER 

DIGITAL-OOMPUTER  SOLUTIONS 

To  solve  actual  numerical  cases ^  the  general  one- 
compartment  (one— phase)  chemical-equilibrium  code  developed 
at  RAND  for  the  method  described  in  Ref.  5  was  extended  to 
cover  the  multicompartment  case.  The  method  is  iterative  and 
starts  with  an  arbltrai'y  composition  of  chemical  species 
satisfying  the  mass— balance  constraints.  By  means  of  the 
assumed  chemical  composition,  a  quadratic  fit  is  made  to 
the  free— energy  function,  which  is  then  minimized.  This 
produces  a  new  chemical  composition  that  also  satisfies  the 
constraints  but  has  a  lower  free  energy.  Tlie  process  is  then 
Iterated.  The  method  is  extremely  efficient  and  for  this 
reason  has  come  i-nto  wide  use  for  solving  chemical— equilibrium 
problems . 

Model  I  (v/hich  has  11  mass-balance  equations  and  50 
chemical  species)  can  be  computed  on  the  IBM  704  to  four  sig¬ 
nificant  figures  in  approximately  one  minute  with  less  than 
20  Iterations.  The  cost  of  computing  the  output  composition 
for  a  given  set  of  Inputs  from  venous  blood  and  air  composi¬ 
tion  is  about  $5.00.  The  currently  used  digital  code  is 
capable  of  handling  systems  consisting  of  10  times  the  number 
of  equations  and  variables.  Hov/ever,  systems  of  such  size  may 
require  a  prohibitive  number  of  iterations  before  an  accurate 
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solution  can  be  achieved.  On  the  other  hand,  systems  of  two 
to  three  times  the  present  size  are  believed  to  be  economi¬ 
cally  feasible.  So  far  as  the  authors  know,  the  solution 
carried  out  to  validate  the  model  (see  Sec.  Vl)  represents 
the  largest  simultaneous,  multiphase  chemical  system  ever 
solved  for  equilibrium. 

For  Model  I,  the  code  used  only  5,000  of  the  52,000  words 
of  high-speed  memory  space  available  (including  space  for 
about  4,000  words  for  the  code  Instructions  and  for  about 
1,000  words  for  the  problem  data,  namely,  the  Input  values  b^^, 
the  energy  constants  Cj,  etc).  In  addition,  the  machine 
has  magnetic— tape  and  drum  memories.  Thus  It  Is  feasible  to 
consider  having  several  human— subsystem  models  within  the 
machine  at  the  same  time;  these  can  be  cyclically  solved  so 
as  to  feed  Input  Information  to  each  other.  For  example,  a 
-respiratory  model  at  the  tissue  level  could  Interact 
with  the  present  external  respiratory  model.  (A  simple  blood- 
flow  simulation  v/ould  be  needed  to  Interconnect  the  two  sub¬ 
systems.)  This  two— subsystem  model  should  give  a  fairly 
complete  simulation  of  the  full  respiratory  function  of  the 
human  body.  One  Immediate  application  might  be  to  the  study 
of  body  reactions  when  breathing  In  confined  atmospheres 
such  as  are  envisioned  for  space  travel. 
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ANALOa-COMPUTER  SIIWLATION 

In  addition  to  obtaining  accurate  numerical  values  using 
a  digital  computer,  the  authors  explored  the  Idea  of  simulating 
the  respiratory  function  (so  far.  Just  the  chemical— equilibrium 
part)  on  the  electronic  circuits  of  an  analog  computer.  For 
this  purpose,  the  RAND  analog  computer  (iRAC)  was  used.  This 
Is  a  differential  analyzer  originally  made  by  the  Reeves 
Instrument  Company,  but  with  many  modifications.  The  various 
elements  of  the  machine  Integrators  (which  are  potentiometers) 
were  Interconnected  by  means  of  a  wired  plug  board  that  con¬ 
stitutes  the  Model  I  "program"  for  the  machine  to  execute.  A 
series  of  dials  can  be  turned  to  vary  the  Inputs  to  the  de¬ 
sired  amounts,  and  the  circuits  of  the  machine  can  be  scanned 
to  read  the  values  of  any  desired  output  species. 

Computer  solutions  of  Model  I  were  undertaken  with  two 
goals  In  mind.  The  first  was  the  obvious  necessity  of  having 
fast  accurate  checks  on  the  mathematical  model;  the  second  was 
the  desirability  of  obtaining.  In  physical  hardware,  a  simu¬ 
lation  of  the  time-dependent  system.  Although  these  two  goals 
are  not  Independent,  the  digital  computer,  with  high  accuracy 
and  reasonably  fast  solution  time,  was  used  to  compute  the 
equilibrium  conditions  and  to  check  the  mathematical  model. 

The  analog  computer,  v/lth  less  accuracy,  but  much  faster  solu¬ 
tion  time,  can  be  used  eventually.  It  Is  believed,  to  simulate 
the  dynamic  aspects  of  the  model.  The  Initial  effort  to  this 
end  v/as  accomplished  In  the  following  way. 
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The  analog  computer,  as  its  name  Implies,  Is  singularly 
well  designed  for  the  simulation  of  complex  time-dependent 
systems.  With  the  advent  of  mathematical  programming  methods. 
It  became  possible  to  simulate  complex  systems  for  which  the 
basic  problem  can  be  stated  In  tenns  of  the  extremlzatlon  of 
a  given  vector  —  for  example,  the  chemical— equilibrium 
problem  of  Model  I  as  stated  above.  By  means  of  the  methods 
of  mathematical  programming,  an  exact  model  of  the  equations 
for  Model  I  was  constructed  on  the  analog  machine.  An  Im¬ 
portant  characteristic  of  this  simulation  Is  that  the  entire 
system  Is  Interconnected  and  hence  Interdependent  In  Just  the 
way  the  hypotheses  of  the  physiological  system  are  Interde¬ 
pendent.  This  means  that  a  change  In  level  of  concentration 
or  activity  of  any  chemical  species  or  component  Is 
Immediately  reflected  throughout  the  system. 

All  the  Important  parameters  of  the  system  were  made 
available  on  the  machine  for  manual  or  automatic  adjustment. 
Once  the  model  had  been  constructed  on  the  machine,  the  Input 
or  known  parameters  were  adjusted  by  means  of  potentiometers. 
The  machine  quickly  computed  the  equlllbrlrmi  values  for  the 
outputs  of  the  chemical  model.  That  Is,  the  model  on  the 
machine,  with  circulating  currents,  came  to  an  equilibrium  In 
which  the  current  level  and  voltage  were  constant,  a 
simulation  of  the  chemical  system.  The  outputs,  or  equilibrium 
values,  could  then  be  read  on  voltmeters  or  appropriate  output 
devices  and  the  comparison  made  to  the  physiological 
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hypothesis.  For  example,  the  effect  on  the  amount  of  water 
in  the  red  blood  cells  was  studied  by  varying  the  amount  of 
salt  in  the  plasma.  This  was  done  by  continuously  changing 
the  setting  of  the  Na"*"  and  Cl~  Input  dials.  Similarly  the 
effect  of  acidity  on  oxygen  pickup  v/as  studied  by  varying 
the  H'*’  Input,  etc.  It  is  quite  possible  that  such  a  simu¬ 
lator  vfith  more  subsystems  would  be  an  excellent  training 
device  to  give  a  medical  student  at  the  controls  a  feel  for 
the  effect  of  this  or  that  change  on  the  functioning  of  the 
human  system. 

To  make  such  a  system  time-dependent  in  the  most 
elementary  way,  the  input  parameters  (potentiometers)  —  which, 
for  example,  set  the  CO2  concentration  in  the  lung  alveolar 
sacs  —  were  made  to  vary  in  a  sinusoidal  fashion  with  a 
frequency  comparable  to  the  breathing  rate.  The  entire  system 
oscillated  according  to  the  interdependent  relations  at  the 
same  rate.  These  oscillations  may  not  correspond  to  nature, 
of  course,  because  a  more  complete  representation  may  require, 
for  example,  the  simulation  of  the  time  lag  encountered  across 
physiological  boundaries,  the  interconnection  of  similar  simu¬ 
lations  of  other  organs,  and  a  simulation  of  the  respiratory 
function  at  cellular  level. 

The  simulation  of  Model  I  required  approximately  50 
computing  amplifiers  and  50  potentiometers.  The  small  analog 
computer  at  The  RAND  Corporation  can  be  made  to  simulate  a 
system  about  tv/lce  the  size  of  the  present  Model  I  —  e.g.. 
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the  full  respiratory  function.  For  larger  systems,  there  are 
machines  available  outside  RAND.  It  Is  also  possible  that, 
as  subsystems  become  fully  validated,  special-purpose  analog 
computers  could  be  built  for  the  sole  purpose  of  simulating 
the  particular  subsystems. 
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VI.  VALIDATION  OF  MODEL  I 

Where  possible^  the  respiratory  model  used  equilibrium 
constants  derived  from  physical  chemical  data  rather  than 
those  derived  from  laboratory  measurements  on  the  biological 
medium  Itself.  The  rationale  for  this  choice  Is  that  data  of 
the  second  type  may  reflect  the  presence  of  other  substances 
and  reactions  and  may  be  valid  only  for  the  normal  or  measured 
conditions.  On  the  other  hand,  the  use  of  basic  data  of  the 
first  type  should  allow  explorations  with  the  model  of  con¬ 
ditions  that .are  far  from  normal  In  respect  to  temperature, 
pressure,  blood  and  air  composition,  etc.,  because  the  more 
basic  equilibrium  constants  can  be  expressed  functionally  In 
terms  of  these  conditions.  The  Important  first  question  to  be 
asked  Is  how  well  a  model  constructed  In  this  manner  conforms 
to  the  actual  physiological  phenomena. 

To  validate  the  model,  the  problem  chosen  was  to  find  the 
equilibrium  amounts  of  each  of  the  species  for  the  Inputs 
shown  In  Table  2  for  an  average  resting  adult  male  (at  1 
atmosphere  and  57^C) .  This  was  accomplished  by  minimizing  the 
sum  of  the  Individual  partlal-molal  free-energy  functions  for 
the  species,  subject  to  the  linear  constraints  for  the  con¬ 
servation  of  mass  and  of  charge,  using  the  electronic  computer 

code  described  In  Sec.  V.  Table  3  lists  the  calculated  values 

(q) 

from  the  model  and  the  observed  values  reported'^'  for  average 
resting  Individuals.  It  can  be  seen  that  the  correspondence 
of  the  calculated  and  observed  values  Is  quite  close. 
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Table  2 

INPUTS  FCU  THE  RESPIRATION  SYSTEM  MODEL® 


''Element'' 

From  Air  (moles 
per  0.051  moles 
air) 

From  Mixed— venous 
Blood  (moles  per 
liter) 

Total 

^2 

a^  =  6.4872  X  lO"”^ 

Vj^  =  6.85  X  10~^ 

b^  =  1.53172  X  10“^ 

CO2 

a^  =  9.27  X  10~^ 

V2  =  2.27  X  10“^ 

bg  =  2.270927  X  10~‘ 

^2 

a^  =  2.44176  X 10“^ 

v^  =  4.57  X  10“^ 

b^  =  2.48546  X  10~^ 

h'*' 

a^  =  0 

Vj^  «»  46.7 

bj^  =  46.70 

0H“ 

a^  =*  0 

v^  -  46.71975 

b^  =  46.71975 

Cl“ 

ae  =  0 

vg  =■  0.0814 

bg  =  0.0814 

Na'*' 

a^  =  0 

v.^  =  O.O8O92 

b^  =  O.O8O92 

as  =  0 

Vg  =  0 . 050 

bg  =  0.050 

HB1“^ 

ag  =  0 

Vg  =  9.09  X  10“^ 

bg  =  9.09  X  10“^ 

®10  “  0 

Vio  “  8.80  X  10~^ 

bio  " 

HP,-" 

®11  =  0 

v^^  -  1.19  X  10“^ 

bn  =  1-19  X  10-^ 

Charge  z 

®12  “  ° 

V12  0 

bl2  ■  0 

®These  data  have  been  taken  from  Ref.  9;  they  represent 
averages  for  resting  adult  males. 

'^The  symbol  B1  denotes  1/4  reduced  hemoglobin  molecule, 
mol.  v;t.  16,500;  symbols  Pp  and  P^  represent  miscellaneous 

protein  in  plasma  and  red  blood  cells,  respectively. 
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Value  from  Model  Observed  Value 
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Some  interesting  features  of  this  equilibrium  model  may 
be  observed  by  following  some  of  the  reactions  and  solutions 
in  Tables  1,  2, and  5*  For  example,  it  is  Important  to  ob¬ 
serve  the  respiration  quotient  (RQ) ,  the  moles  of  COg  out 
divided  by. the  moles  of  Og  in.  The  moles  of  COg  out  may  be 
obtained  by  subtracting  the  value  for  COg  input  (Table  2) 
from  0.0326879  times  the  Xg/  x  value  for  COg  output  (Compartment 
I,  Table  3)*  If  a  similar  subtraction  is  carried  out  for  Og, 
the  results  give 


„  o^ooitiM  ,  0.86, 
0. 002067  "l 


a  value  well  within  the  range  for  resting  individuals. 
Similarly,  the  per  cent  of  hemoglobin  saturation,  as  de¬ 
termined  by  the  model,  is  the  reasonable— appearing  value 


100  Xgg 

X28  +  Xgg 


0.878 

■o':'oo'905 


96.6  . 


Inar much  as  the  carbamlno  reaction  involving  the 
combination  of  hemoglobin  and  COg  v/as  not  represented  in 
Model  I,  it  was  necessary  to  compensate  for  this  lack  of 
a  COg  sink  in  arterial  blood  by  reducing  the  total  venous  COg 
input  by  a  comparable  amount.  Consequently  the  Vg  inprut  value 

_3 

indicated  is  1  *  10  moles  less  than  the  average  true  value 
for  total  COg  in  venous  blood  to  compensate  for  the  lack  of 
carbamlno  representation.  Even  though  the  carbamlno  reaction 
Involves  a  small  part  (approximately  5  pei*  cent)  of  the  COg 
carried  by  the  blood,  it  accounts  for  a  much  higher  portion 
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I 

(approximately  50  per  cent)  of  the  CO2  ti’nnnferred  by  the 
blood  during  the  respiratory  cycle-  A  model  of  the  next 
degree  of  sophistication  must  therefore  account  for  this  im¬ 
portant  reaction.*  The  representation  of  the  amounts  of  the 
other  gaseous  constituents  In  Compartment  I  —  njunely,  N2  and 
HgO  —  appears  to  be  nicely  In  line  with  the  Inputs  and  with 
reported  values. 

In  the  cases  of  columns  x^,  x^,  x^,  which 

show  the  amounts  of  0^,  COg,  and  N2  dissolved  in  plasma  and 
red  cells,  respectively,  it  will  be  noted  that  there  are 
some  rather  wide  variations  between  the  values  calculated 
by  the  machine  and  the  observed  values.  This  is  most 
apparent  for  ©2  and  Ng  dissolved  in  red  cells.  Wiere  is 
a  possibility  that  the  model  results  may  be  more  nearly 
correct  than  the  observed  values  in  these  cases.  Such  very 
small  quantities  of  dissolved  gases  are  difficult  to  de¬ 
termine  in  the  laboratory.  Then,  too,  it  is  possible  that 
small  amounts  of  hemoglobin  remained  in  the  analytical 
samples,  which  combined  with  oxygen  and  thereby  made  the 
solubility  appear  to  be  greater  than  it  is  in  the  complete 
absence  of  hemoglobin.  Although  the  number  of  moles  of  these 
dissolved  gases  Is  different  for  plasma  and  red  cells,  their 
mole  fractions  are  identical  in  both  comjjartments  in  the 
machine  solution.  This  is  so  for  all  species  appearing  except 

*For  an  e::oc*llent  description  of  the  physiological  role 
of  carbamlno  compounds  in  blood,  see  Ref.  5,  pp.  87—88. 
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for  those  charged  species  that  are  influenced  by  the  Oibbs— 
Donnan  equilibrltm.  It  should  be  noted  that  the  solubility 
coefficients  used  in  the  model  are  for  water>  not  for  plasma 
and  red  cells. 

It  is  also  interesting  to  note  that  the  model  creates 
water  in  the  three  compartments  from  the  Inputs  and  Olf* 
operating  through  the  reaction 

+  OH"  HgO  . 

The  amount  of  water  formed  is  Influenced  slightly  by  the 
appropriate  demands  and  supplies  of  and  OlT"  by  the  com¬ 
peting  reactions. 

The  greatest  difficulties  in  connection  with  the  valida¬ 
tion  of  the  model  arose  in  connection  with  the  representation 
of  acidity  and  concentrations  of  and  colximns  Xq, 

XgQ#  and  in  plasma  and  red  cells.  Acidity  of  blood 

samples  is  usually  determined  by  the  use  of  a  pH  meter  and  a 
glass  or  similar  electrode.  It  is  common  to  see  the  ex¬ 
pression 

pH  E  -loSjo  =  “^°®10  ®H‘*‘  * 

However,  in  recent  years  conclusive  arguments  have  been  pre¬ 
sented,  based  on  both  theoretical  and  practical  electrolytic- 
cell  considerations,  that  pH  values  are  not  measures  of 
hydrogen— ion  concentration  or  activity.  Rather,  pH  determina¬ 
tions  are  a  reliable  measure  of  the  "acidity"  of  the.unknovm 
solution  in  relation  to  an  empirical  standard  solution  having 
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the  same  general  characteristics.  Including  temperature. 
Conventional  pH  measurements  are  of  use  primarily  for  control 
and  gross  comparison  purposes,  but  are  not  suitable  for 
thermodynamic  calculations  or  mass— action  law  use.  There  are 
not  even  really  satisfactory  methods  for  approximating 
or  [Old  from  pH  values  or  vice  versa. 

The  problem  this  presents  In  the  validation  of  the 
operation  of  the  model  la  to  determine  If  the  model  calcu¬ 
lates  [hT  and  [on  accurately  when  the  only  laboratory 
measure  of  these  concentrations  for  comparison  are  pH  values 
for  plasma  and  red  cells.  One  approximation  suggested  In  the 
literature la  that 

pH  measured  i  ApH  =  • 

when  [hT  la  determined  from 

-  [Hi  [0H-]  ,  . 

[HjO] 

where  la  established  by  other  means,  such  as  conductivity 
measurements.  It  will  be  noted  that  If  the  reported  pH 
values  for  plasma  and  red  cells  are  taken  equal  to  — log^^  [hT, 
and  — log^Q  Is  determined  by  subtraction  from  log 

(which  la  —15*55  at  57°C),  then  the  resulting  concentration 
product  does  not  equal  K^;  this  Indicates  that.  In  blood, 

pH  measured  /  [hT. 

However,  If  one  assumes  that  the  acid— base  reactions  finally 
Involve  only  the  vrater  contained  In  the  plasma  and  red— cell 
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material,  and  If  one  reduces  the  given  pH  values  (and 
calculated  [on  values)  to  correspond  to  those  for  which 
water  only  is  present,  then  the  values  so  obtained  for  real 
blood  satisfy  the  several  mass— action  relationships  quite 
nicely  and  compare  well  with  the  values  calculated  by  the 
model  and  listed  In  Table  2.  This  correction  was  applied  to 
the  values  of  and  [oH~3  listed  as  observed  values  In 

Table  2  and  as  approximated  from  the  reported  pH  measurements. 

The  treatment  of  species  and  x^q,  sodium 

Ion  and  proteins  In  plasma  and  potassltim  Ion  and  proteins  In 
red  cells.  Is  described  In  Appendix  B.  Species  and  Xg^# 
HCOj  in  plasma  and  in  red  cells.  Is,  of  course,  a  cinjclal  Ion 
in  the  blood  both  for  COg  transport  and  for  buffering.  It  Is 
created  by  the  model  In  Compartments  II  and  III  by  the  reaction 


COg  +  OH"  HCOj 


In  the  face  of  the  demands  and  supplies  of  these  Inputs  from 
competing  reactions.  Bicarbonate  Ion  Is  said  to  hydrolyze 
(see  p.  219  of  Ref.  ll) l 

HCOj  +  HgO  HgCOj  +  Otf" ; 

and  it  can  ionize  further i 


HCO,  i±  C0“  +  H"^  . 
P  p 


The  combination  of  these  two  reactions  gives 


2HCO”  +  HgO  HgCOj  +  C0“  +  H"^  +  OH"  . 
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Species  Xgg,  Xgy  —  namely  H2C0^  and  CO^  in  plasma 

and  red  cells  —  with  their  appropriate  thermodynamic  func¬ 
tions,  were  introduced  in  the  model  to  explore  their  effects 
on  concentrations  of  HCO^,  and  Off”.  The  effects  of  HCOj 

hydrolysis  and  second— stage  ionization  on  [HCO^  and  acidity 
appear  to  be  minor. 

The  concentrations  of  species  H2C0^  and  CO*  are  not  re¬ 
ported  in  the  literature.  These  species  viere  retained  in 
the  model  to  Illustrate  the  possibilities  of  using  the  model 
for  explorations  of  this  sort.  In  addition,  it  has  been 
noted  that  the  model  is  extremely  and  correctly  sensitive  and 
responsive  to  changes  in  the  inputs  or  in  the  environment 
(e.g.,  temperature  and  pressure). 
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VII.  POTENTIALITIES  AND  LIMITATIONS  OF  MODEL  X 

The  authors  hasten  to  admit  that  the  model  described  is 
deficient  in  many  respects,  for  it  represents  but  a  first 
attempt  to  determine  the  feasibility  of  programming  a 
complicated  physiological  system  in  this  manner.  Because  the 
Interest  of  this  particular  study  is  centered  on  the  equili¬ 
brium  phenomena  of  the  system's  chemical  aspects  rather  than  on 
the  physical  flow,  mixing,  etc.,  sevaral  expedients  have  been 
used  to  av  .id  incorporating  a  representation  of  the  physical 
and  transient  aspects  of  the  entire  external  respiratory 
system  in  the  model.  For  example,  the  mechanism  of  breathing 
with  its  complex  intermittent  flow,  mixing,  and  diffusion 
effects  has  been  avoided  for  this  initial  study  in  several 
different  ways: 

1.  By  making  the  set  of  input  "elements,”  a^,  ...,  a^, 
representative  of  a  large  excess  volvime  of  air  having  the 
observed  composition  of  the  gases  of  the  pulmonary  alveoli  of 
a  normal  resting  male  at  sea  level.  The  assumption  introduced 
by  this  procedure  is  that  the  concentration  of  gases  in  the 
alveoli  of  resting  individuals  does  not  change  significantly 
during  the  breathing  cycle. 

2.  By  relating  the  required  Incremental  input  of  "fresh" 
air  of  atmospheric  composition  to  the  Increments  of  gaseous 
"elements"  needed  to  convert  a  unit  volume  of  mixed— venous 
blood  to  arterial  blood.  The  "fresh  air"  Inputs,  a^,  ...,  a^. 
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per  liter  of  blood  may  be  estimated  on  this  basis  when  this 
Information  Is  not  provided  by  the  source  data. 

5.  By  programming  the  computer  to  sample  the  concentra¬ 
tions  of  ©2  and  CO^  In  either  the  arterial  blood  or  sac  air 
compartments  and  calculate  by  a  series  of  Iterations  the 
amounts  of  "fresh"  air  of  given  composition  required  to  make 
these  lung  or  blood  values  correspond  to  previously  established 
standards. 

All  three  of  these  techniques  have  been  used  and  are 
satisfactory  within  their  limitations.  Method  (2)  above  was 
used  In  obtaining  the  values  shown  In  Tables  2  and  5.  These 
devices  avoid  the  necessity,  at  this  early  stage,  of  combining 
a  physical  model  of  breathing  with  the  present  chemical  model, 

which  may  require  the  Introduction  of  a  number  of  concepts  In 

» 

addition  to  those  needed  for  the  chemical  model. 

Another  assumption  Implicit  in  the  model  Is  that  the 
arterial  blood  leaving  the  lungs  comes  to  equilibrium,  through 
the  alveolar  and  capillary  walls,  with  the  gases  In  the  alveoli. 
As  noted  earlier,  realization  of  this  equilibrium  situation 
appears  to  be  closely  approximated  by  resting  Individuals. 

Hence,  for  the  one-cycle  case  used  to  validate  the  model 
(approximately  0.55  liter  of  arterial  plasma  plus  approximately 
0.45  liter  of  arterial  red  cells  In  equilibrium  with  each  other 
and  with  approximately  0.033  moles  of  alveolar  gases,  all 

« 

One  mathematical  representation  of  the  breathing 
mechanism  can  be  found  In  Ref.  12. 
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subject  to  the  mass-balance  constraint  provided  by  the  "elements" 
of  one  liter  of  mixed-venous  blood  plus  G.O31  moles  of  alveolar 
gases),  the  assumption  is  made  that  equilibrium  is  reached 
between  the  solutions  within  the  red  cells  and  the  plasma.  A 
comparison  of  the  results  obtained  from  the  model  with  those 
obtained  from  laboratory  data  on  blood  samples  shows  that  this 
assumption  does  not  appear  to  introduce  any  serious  error  for 
resting  individuals. 

The  multistaged  ionization  of  hemoglobin  and  the  other 
protein— like  blood  substances  is  not  represented,  and  the 
Important  carbamlno  reaction  is  only  implicitly  represented, 
in  Model  I.  In  Appendix  D,  an  expanded  system  (Model  II)  is 
shown,  demonstrating  how  these  factors  may  be  incorporated. 

At  the  present  stage,  the  model  is  one  of  instantaneous 
reaction,  with  time  not  explicitly  Introduced  as  an  independent 
variable.  This  m.eans,  of  course,  that  the  present  representa¬ 
tion  can  deal  only  with  features  of  the  respiratory  system 
that  are  capable  of  treatment  by  reversible  thermodynamics; 
however,  this  limitation  still  allows  many  matters  of 
importance  and  Interest  to  be  explored.  In  general,  if  the 
number  of  relatively  slow  chemical  reactions  entering  the  model 
is  not  large,  it  may  be  possible  to  develop  a  dynamic  model  by 
introducing  some  simple  delays  in  the  electrical  analog  of  the 
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« 

current  model. 

It  will  also  be  necessary  to  develop  and  Introduce  the 
feed-back  loops  by  which  changes  In  COg  concentration,  and 
possibly  hydrogen-ion  concentration.  Influence  the  ventilation 
rate  and  the  rate  of  blood  flow  through  actions  on  the 
cardiovascular  system.  Consideration  must  also  be  given  to 
the  bo’indary  conditions  essential  for  the  continued  functioning 
of  the  subsystem.  One  Important  limit  Is  the  maximum  osmotlo 
load  tolerated  by  cells.  Fortvinately,  the  model  Is  nicely 
suited  to  the  Incorporation  of  this  limit. 

In  spite  of  these  difficulties.  It  should  be  possible, 
through  the  use  of  models  of  this  type,  to  demonstrate  the 
detailed  Interactions  of  the  human  respiratory  system  as  the 
external  environment  changes  and  as  the  Inputs  from  related 
subsystems  change  (e.g..  In  metabolism).  Changes  In  total 
respiratory  gas  pressure  or  relative  mixtures  will  be  reflected 
In  the  model  by  the  appropriate  changes  In  oxyhemoglobin 
concentration,  alkali  reserve,  blood— flow  rate,  breathing  rate, 
plasma  composition,  etc.  Because  the  energetics  are  an 
Integral  part  of  the  model,  changes  In  work  load,  temperature. 


# 

In  highly  dynamic  situations,  such  as  when  the  respira¬ 
tory  system  Is  stressed  by  exercise,  there  Is  good  reason  to 
believe  that  the  controlling  reactions  will  be  those  Involving 
membrane  transfers.  There  are  three  such  transfers  In  the 
external  respiratory  system.  Involving  (l)  the  alveolar  tissue 
In  the  lungs,  (2)  the  red-cell  membrane,  and  (5)  the  body— cell 
membrane.  The  Introduction  of  the  appropriate  kinetic  equations 
for  any  or  all  of  these  three  barriers  may  suffice  to  explain 
all  the  chemical  time  effects.  The  remaining  reactions  occur  as 
homogeneous  reactions  In  dilute  solution  or  In  the  gas  phase  and 
are  thus  generally  very  fast  compared  to  membrane— diffusion 
processes. 
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and  pressure  Imposed  on  the  human  system  should  be  reflected 
by  the  appropriate  and  detailed  response  within  the  modeled 
respiratory  system.  Changes  In  fields,  either  gravitational, 
electromagnetic,  or  electrostatic,  can  also  be  Introduced 
through  the  thermodynamics.  Abnormal  changes  In  any  of  the 
constituents  of  the  respiratory  system,  as  brought  about  by 
pathological  conditions  or  by  unusual  environments,  should 
elicit  the  appropriate  response  In  the  other  constituents. 
Including  the  energetic  aspects.  Thus  an  external  respiratory 
submodel  after  further  development  should  prove  most  useful 
for  diagnosis  and  for  predicting  the  possible  effects  of 
vinusual  stresses  on  the  respiratory  system  without  endangering 
human  life. 


RM-2519 

57 


Appendix  A 


A  REVIK'/  OF  CERTAIN  CHEMICAL  THERI'IODYNAMIC  CONCEPTS 

For  the  benefit  of  readers  who  are  not  familiar  with 
chemistry,  we  shall  review  some  of  the  basic  concepts  used 
In  our  study,  especially  the  relation  of  the  mass--actlon 
law  to  minimum  free  energy.  Suppose  a  compartment  contains 
various  species  such  as  HgO,  OK”,  COg,  HCO^,  when 
COg  Is  dissolved  In  the  solution.  The  amounts  of 
these  various  constituents  depend  on  the  total  amounts  of 
water  and  carbon  dioxide  originally  "dumped"  Into  the  compart¬ 
ment,  but  their  relative  concentrations  to  each  other  satisfy 
certain  conditions  that  are  Independent  of  the  Input  amounts. 

In  general,  the  concentrations  (or  activities)  of  the  species 
at  equilibrium  have  the  property  that,  independent  of  the  In¬ 
put  quantities,  one  or  more  ratios  of  products  of  concentra¬ 
tion'  among  them  have  fixed  values.  These  nonlinear  conditions 
are  called  mass— action  relationships. 

For  example.  If  we  let  x  denote  the  total  number  of 
molecules  In  the  mixture,  so  that 


^  “  ^HgO 


^0H~  ^COg  ^HCOj  * 


where 


-  the  number  of  un— Ionized  water  molecules. 


x^„—  -  the  nimiber  of  Ions  of  Orf”,  etc,. 


and  define 


AUrjr-C  'jxy 
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to  bo  the  ooncontratlons  of  these  species  In  the  mixture,  then 
the  niasB--eotl6n  law  states  that 


^OS-)  (v)  ^ 


whojTe  Kj  has  a  fixed  value  (at  constant  temperature  and 
pressure),  called  the  equlllbrltun  constant  for  the  reaction 


H'*'  +  Oir  ^  HgO  . 


Even  If  no  COg  were  present,  or  If  a  salt  were  added  to  the 
nlxtuM,  this  same  relationship  would  hold  despite  the  fact 
that  the  Individual  concentrations  forming  the  above  ratios 
could  be  changing  drastically.  Another  mass-action  relation¬ 
ship  that  holds  among  the  concentrations.  Independent  of  the 
amounts  of  Inputs,  Is 


where  Kg  Is  the  equilibrium  constant  for  the  reaction 

COg  +  HgO  ^  HCO“  +  . 


^  HCOJ  +  H"*"  . 
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Ihe  oonvantlonal  method  for  determining  the  equilibrium 
oompooltlon  Is  to  find  the  values  of  the  five  unkno^ms  Xj  that 
satisfy  five  equations  —  two  of  which  are  nonlinear  and  are 
the  two  mass— action  conditions  given  above,  while  the  remain¬ 
ing  three  are  linear  and  express  the  mass  balance  In  relation 
to  the  original  Inputs,  for  example. 


^OH  *HCOj 


Input  Olf, 
Input  H^, 


^COg  *HCO“ 


Input  COg  • 


As  noted  earlier  In  this  research  memorandum,  we  have  found 
It  more  convenient  to  bypass  the  mass-action  relationships  by 
finding  the  equilibrium  composition  through  direct  minimiza¬ 
tion  of  the  free— energy  function^  subject  only  to  the  mass- 
balance  conditions. 

In  general,  for  a  chemical  reaction  represented  by 


+  Ag  Pg  +  •  «  «  P^  +  Xg  Pj^  +  •  •  • ,  (T) 

the  corresponding  mass-action  relationship  Is 


■■■ 


t 

where  the  Pj  are  the  reactant  species;  the  Pj  are  the  product 

I 

species,  the  Xj  and  Xj  are  the  corresponding  niunbers  of  moles 
of  each  species,-  K  Is  the  equlllbrlvun  constant;  and  the 
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Nj  and  Nj  are  the  concentrations  of  the  constituents  Pj  and 
pj.  We  may  also  state  the  mass-action  law  In  logarithmic  form 
by  taking  logarithms  of  both  sides;  thus. 

In  K  - 


t  t  I  - 

In  Nj  +7^2  In  Ng  +  ... 


lx  In  +  Xg  In  Ng  + 


..]  . 


(9) 


It  Is  Important  to  note  that  a  chemical— reaction  equation# 
denoted  by  (7)^18  simply  f  relation  among  the  colximn  vectors 
Pj  of  Input  elements  associated  with  each  species  Pj, 

Thus  If  we  let 


we  have#  corresponding  to  (7), 


X^  P^  +  Xg  Pg  +  ...  -  X^*  P^'  +  Xg’  Pg»  +  .  .  .  .  (lO) 

This  Is  easily  seen  from  an  example.  Suppose  H"*"  and  Orf”  are 
considered  Input  "elements"  forming  species  OlT*#  and  HgO. 
We  form  the  tableau 


Input 

Species 

OH 

HgO 

Element 

1 

0 

1 

oir 

0 

1 

1 

where  the  entry  a^^j  In  row  1  and  column  j  Is  the  quantity  of 
the  Input  element  1  found  In  the  species  J.  It  Is  clear  that 
the  chemical  reaction 


+  Oir  ^  HgO 


i 
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Is  simply  a  statem^crnt  about  the  correspondlnc  column  vectors 
in  the  tableau,  namely 


’l' 

'o' 

’l" 

+ 

as 

0 

1 

_r 

We  aliall  novr  shov;  that  the  mass--action  law  holds  for  the 


values  of  Xj  that  minimize 


n 


F(X)  -  Z  X  P. 

v^ere  In  particular  Pj  Is  given  by 


F.  =  RT  (c  +  In  )  . 
J  .)  X 


(11) 


(12a) 


As  noted  earlier,  this  expression  for  Fj  Is  essentially  an 
approximation  for  dilute  solutions  and  Ideal  gases  of  the  more 
general  expression  (l2b).  It  will  be  assumed  unless  otherwise 
stated  that 


Fj  =  RT  (Cj  +  In  aj)  , 


(12b) 


where  each  aj,  the  "activity,"  Is  some  function  of  the  con¬ 
centration  Xj/  x;  that  Is  to  say,  a^  (and  hence  Pj,  the  free 
energy  per  mole  of  the  J— th  constituent)  remains  unchanged  If 
all  quantities  Xj  are  increased  proportionally.  V/e  now 
establish  the  following  purely  mathematical  theorem: 

Theorem :  If  for  every  J ,  Fj  =  Pj(Xj/  x)  Is  a  function 

of  Xj/  X  only,  then 

^F 
■j 


sr:  “  "j' 


J  —  1,  2,  •••,  n* 


(13) 


1 
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Because  of  {\J>) ,  Fj  Is  commonly  called  the  partial  molar 
free  enerey  of  the  .1— th  constituent.  To  prove  this  relation, 
note  first  that  the  free  energy  F(x)  Is  a  homogeneous  function 
of  degree  one  In  the  variables  (x^,  Xg,  •••, 

homogeneous  functions  of  degree  one  It  Is  generally  true  that 


F(X)  =  Z  X.  ^  .  (14) 

J=1  ^ 

This  relationship  (l4)  Is  known  as  Euler's  theorem  for  homo¬ 
geneous  forms  of  the  first  degree.  In  words,  it  says  that 
total  free  energy  can  be  determined  ^  If  3F/^Xj  Is  the 
contribution  of  the  J— th  constituent  per  mole. 

The  foregoing  theorem  Is  easy  to  prove;  we  shall  establish 
It,  as  follows:  From  (l2b),  the  assiamed  homogeneity  of  degree 
zero  for  Fj  for  all  t,  we  have 

tF(x^,  •••»  =  F(tx^,  tx2,  .••,  4Xj^)  . 

Hence,  taking  the  partial  derivative  of  both  sides  with  re¬ 
spect  to  t,  we  obtain 


n^F(tx-,  tx^, 

.)  =  Z  - 


d(tXj) 


Setting  t  «  1  yields  (l4). 

To  establish  (15),  we  now  assume  that  F(x)  Is  given  by 
(11),  where  Fj  -  Fj(xj/  x) ,  and  take  the  partial  derivative  of  F(X) 
with  respect  to  Xy  This  yields 
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(15) 


where  we  have  replaced  bytP^/dx.  since  Ox/^Xj^-  1. 

We  now  show  that  the  suminatlon  term  vanishes.  Multiplying 


both  sides  by  Xj  and  summing  on  J  yields,  by  relations  (l4) 


and  (11), 


Since 


n 


X.  5  0,  and  S  x.  >  0  , 
^  J-1  J 


the  second  factor  of  the  product  must  vanish  and  (15)  follows 


from  (15). 


liavlng  established  that  Pj  a  dP/^j,  let  us  turn  to  the 


main  problem  of  minimizing  the  free— energy  function  F(X) 
subject  to  mass-balance  constraints. 


n 


(16) 


•  •  •  ^  171  • 


We  now  revlev  the  Lagrange  method  for  minimizing  a  general 
function  P(x)  subject  to  m  general  constraints,  (J),  (x)  -  0* 
and  later  specialize  the  functions  The  first  step  la  to 

assign  unknown  multipliers  to  the  functions  <|*i(X)  and  seek 
the  unconstrained  minimum  of  the  function 


0(X)  -  P(X)  -  Tj  (t)i(x)  -  ...  -  •ir^<I»n,(X)  .  (1?) 
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If  for  some  choice  of  tTj^,  the  point  X  =  X°  where  the  nlnimunj 
of  Cr(X)  occurs  happens  to  satisfy  the  constraints  -  0, 

then  this  unconstrained  minimum  point  X°  for  Q(X)  Is  the  con¬ 
strained  minimum  point  for  P(X).  To  see  this,  note  that 

min  a(X)  -  0(X°)  •  F(X°)  , 

and  that  for  ull  other  X  that  satisfy  (l)j^(X)  0  we  have 

F(X)  2  F(X°)  because 

F(X)  »  a(X)  ^  min  a(X)  -  F(X°)  . 

When  P(X)  possesses  partial  derivatives,  and  assumes  Its  mini¬ 
mum  at  some  Interior  point  of  the  domain  of  definition.  It  can 
be  shown  that  such  multipliers  always  exist.  We  now  specialize 
(Jj^(X)  to  be  linear  as  given  by  (l6).  In  this  case,  we  may 
rewrite  Q(X)  by  collecting  all  terms  with  common  Xj,  obtaining 

a(X)  =  F(X)  -  2(TrP^)  Xj  +  TTb  ,  (l8) 

where 

TT  (1^2 ,  Tr2»  •••» 

Is  a  row  vector, 

®mj^ 

Is  the  coliimn  vector  of  coefficients  associated  with  species 
J  defined  earlier,  and 

b  “  (b^,  bg,  ...,  bj^) 

Is  the  column  vector  of  constants.  The  coefficient  of  Xj  Is 
the  constant  irPj  and  depends  on  the  choice  of  ir.  At  an  un¬ 
constrained  minimum  of  G(X),  X  satisfies 

^  -  (irPj)  -  0 

for  any  choice  of  tt.  Also  for  any  given  chemical  reaction 


I 

I 
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(7)  we  may  multiply  the  corresponding  vector  relation  (10) 
by  v,  obtaining 

>‘l(TrPi)  +  +  ...  “  X^(irP^)  +  (I9) 


vfhlch  holds  for  any  choice  of  tt.  Thus  substituting 
dF/^Xj  "  have 


dx. 


+  "Kr 


6p 

bx.^ 


bx[ 


+  K 


<SL 

bx' 


(20a) 


which  holds  for  the  X  that  minimizes  0(X)  for  any  choice  of 
ir.  Since  F  has  a  special  form  (ll)  that  satisfies 

J  J 

we  also  have,  for  any  such  X, 

+  Xg  F^  t  •  •  •  ^2  ^2  *  *  (^®^) 

If  we  further  set 


Fj  *■  RT  (cj  +  In  Oj)  , 

then  upon  substitution,  rearranging  of  terms,  and  dropping  of 
the  common  factor  RT,  we  get  the  relationship 

I  I  I  I 

^l*'!  ^2^2  '  ^l*'!  ^2*^2  '  ' ' ' 

(20c) 

/•III 

-  (Xj^  In  02  +  X2  In  02  +  ...)  -  (Xj^  In  o^  +  X2  In  Og  +  ...) 

for  all  X  that  minimize  Q(X)  for  any  choice  of  ir.  In  pai>- 
tlcular,  for  that  ir  that  yields  X  =  X°  satisfying  the  mass 
balances  (J)^(X)  »  0  and  thereby  minimizing  the  free  energy 
F(X),  this  same  relationship  must  hold.  But  this  holds  even 
If  vxe  start  out  with  a  different  set  of  b^  values  since  the 
expression  Is  Independent  of  b^.  If  we  define  In  K  to  be 
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In  K  —  ^2^2  *  *  *  ^  ””  ^2^2  *  *  * )  ^ 

then  (20c)  establishes  the  mass-action  laif  In  logarithmic 
form  and  in  terms  of  activities.  If  the  activities  are  given 
by  Sj  *•  Xj/  X,  then  this  is  the  same  as  the  mass— action  law 
stated  earlier  in  equation  (9)»  If 

V  * ' 

where  Yj  Is  some  constant,  this  also  yields  the  mass— action 
law  stated  earlier  with  an  adjusted  value  for  In  K.  In  this 
manner,  the  Cj  values  from  the  free— energy  function  can  be 
used  to  define  the  equilibrium  constants.  The  converse  is 
shown  below. 


DETERMINING  Cj  VALUES  FROM  EQUILIBRIUM  CONSTANTS 

As  was  noted  in  the  text,  the  equilibrium  constants  for 
many  common  reactions  are  tabulated  in  physical  chemical  tables. 

V/e  wish  to  discuss  how  to  adjust  the  Cj  values  so  that 
the  equilibrium  constants  are  directly  applicable  to  the  model; 
to  be  precise,  we  note  that  for  those  x  that  satisfy  the  mass 
balances,  the  expression  for  F(x)  may  be  replaced  by 


n  X  ra 

F(X)  =  Z  X.  (c,  +  In  -J-)  -  Z  k.b.  ,  (22) 

J=1  J  J  X  1=1  1  1 


;-/here  the  are  constants  and  the  Cj  bear  a  direct  relation¬ 
ship  to  the  equilibrium  constants.  The  numbers  k^  were 
selected  so  that  vfhen  multiplied  by  the  1— th  mass— balance 
equation  and  the  m  equations  summed  and  subtracted  from  P(X), 
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m  of  the  new  coefficients  of  the  Xj'S/  denoted 
vanish.  In  place  of  the  original  function  P(x) 


# 


,  would 
the  function 


G(X)  -  Z  x.(c,  +  In  i)  (25) 

J  J  X 

was  minimized  subject  to  the  mass-balance  constraints.  Since 

the  two  functions  differ  by  a  constant, 

m 


the  value  of  X  ■  X°  that  minimizes  P(X)  also  minimizes  G(X). 

Any  m  of  the  c’j  can  be  made  to  vanish  by  suitably 
choosing  the  values,  provided  that  the  square  array  of  co¬ 
efficients  associated  with  these  J  Is  nonsingular,  l.e., 

provided  that  their  determinant  /  0. 

For  Model  I,  the  c’j  values  associated  with  species 
and  Oir"  In  Compartment  II  were  made  eimong  others  to  vanish. 


”  °0H  ” 

Now  the  equilibrium  constant  for  water  does  not  depend  on 
whether  we  use  Cj  or  values  for  Its  definition  In  (21), 

Its  value  being  the  same  since  It  must  yield  the  same 
rlght-4iand  member  of  (20c).  Hence,  for  Compartment  II  si)ecles. 


“  °0H~  “  * 

In  other  words,  the  value  of  0*^2  0  directly 

obtained  from  the  equilibrium  constant  for  water.  In  Appendix 
C,  the  entire  set  of  n  —  m  remaining  Cj  values  Is  evaluated 
In  this  manner  from  equlllbrliam  constants. 


68 


Nc  attempt  waa  made  In  the  construction  of  Model  I  to 


n  O 


determine  the  absolute  Cj  values  aa  defined  In  terms  of 
values.  All  Cj  values  shown  are  the  relative  Oj  v^alues  as 
obtained  from  equlllbrltun  constants. 


■ 


Appendix  B 
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CONSTRUCTION  OP  THE  MODEL  —  INPUTS 


The  theoretical  background  given  In  Appendix  A  makes 
possible  a  more  detailed  description  of  the  construction  and 
operation  of  the  model  shovm  In  Table  1.  Hie  first  block  on 
the  left,  labeled  "Inputs,"  contains  the  "elements"  from  which 
the  molecular  species  In  the  three  compartments  —  "Air  Out," 
"Arterial  Plasma,"  and  "Red  Cells"  —  are  constmicted.  These 
Input  "elements"  have  the  values  In  moles,  b^^,  ...,  that 

are  present  per  unit  of  time  In  mixed— venous  blood,*  plus  the 
values  In  moles  of  fresh— air  gases,  a^,  ...,  a^,  entering  the 
lung  sacs  per  unit  of  time.  For  the  validation  of  the  model, 
one  liter  of  average  mixed— venous  blood  was  used  as  an  Input, 
plus  /X/  0.051  moles  of  degraded  "fresh"  air. 

Inputs  b^,  bg,  b^,  a^,  a^,  a^,  are  obvious.  The  Input 
"elements"  In  rows  4  and  5,  and  0H~*,  provide  a  convenient 
way  to  form  water  In  the  several  compartments  and  serve  as  a 


A  troublesome  data  problem  arose  In  attempting  to 
validate  the  model.  Even  among  normal  Individuals,  the  varia¬ 
bility  of  blood  composition  Is  great,  and  In  addition  It  Is 
necessary  for  validation  to  have  laboratory  data  for  venous 
and  arterial  blood  on  a  completely  comparable  basis. 

Fortunately  for  our  purpose,  the  U.S.  Air  Force  has  prepared 
compilations  of  biological  data  In  this  form.  The  data  for 
Inputs  and  for  comparison  with  outputs  were  obtained  from 
Refs.  9  and  15* 

«  * 

As  explained  elsewhere,  the  air  Inputs  may  also  be 
obtained  by  using  a  large  amount  of  air  of  lung— sac  composition 
or  by  programming  the  computer  to  search  for  the  proper  amount 
of  "fresh"  air  to  meet  blood  or  lung  criteria. 
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mass  source  and  sink  for  the  reactions  Involving  and  Olf”, 
for  example, 

CO.  +  0H“  ^  HC0“  . 

It  will  be  noted  that  HC0“  Is  not  one  of  the  Input  "elements, 
although  It  Is  an  Important  constituent  of  blood.  The  model 
creates  this  species,  however,  through  the  reaction  given 
above.  To  provide  for  this  HCO^  formation  In  terras  of  mass 
and  charge,  the  Input  values  of  H*"  and  Olf*  are  set  different 
by  the  amount  of  HCO^  that  Is  present  In  mixed— venous  blood. 
The  COg  Input  bg,  which  Is  also  Involved  In  production  of 
HCO^,  Is  the  total  COg,  Including  that  present  In  combined 
form.  It  Is  not  directly  Involved,  of  course.  In  the  charge 
balance.  The  mass-difference  effects  In  terms  of  H"*"  and  OH” 
caused  by  the  slight  basicity  of  the  blood  Is  so  small, 
approximately  1  x  10~^®  moles,  that  It  does  not  require 
separate  mass  accounting  at  this  stage. 

"Elements”  Na"*”  and  K"*",  In  rows  7  and  8,  are  Introduced 
at  present  essentially  as  fixed  values  In  plasma  and  red 
cells,  respectively.  Inasmuch  as  entries  for  these  species 
occur  only  and  separately  In  these  two  compartments,  and  are 
not  allowed  to  migrate  In  the  model.  In  real  life,  of  course 
there  Is  some  In  plasma  and  Na"*"  In  red  cells,  "niere  are 
also  other  cations  In  both  compartments,  as  established  by 
the  selective  screening  action  of  the  red— cell  membrane.  If 
and  when  the  data  may  be  available,  this  lon-screenlng 
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mechanism  of  the  membrane  can  be  Incorporated  In  the  model  as 
a  series  of  chemical  equations  with  their  accompanying  thermo¬ 
dynamic  functions.  At  the  present  stage,  however,  all  cations 
(except  H^)  In  plasma  and  red  cells  are  stated  In  terms  of 
equivalent  Na"*"  and  K^,  respectively. 

The  "element"  In  row  S,  designated  by  the  shorthand  HBl"”, 

represents  the  total  hemoglobin  present  In  red  cells.  It  Is 

Introduced  at  this  stage  as  moles  of  one— fourth  hemoglobin 

molecule,  molecular  weight  approximately  16,500  (one  Og# 

one  Fe  per  unlt).^^^  The  present  model  does  not  allow  for 

the  hydrogen  Ionization  of  hemoglobin.  The  aggregating 

^  l41 

techniques  suggested  by  Alberty'  '  provide  a  method  for  In¬ 
corporating  the  multlstaged  Ionization  of  hemoglobin  Into  the 
model  with  the  addition  of  a  modest  number  of  columns. 
Unfortunately,  the  data  provided  for  the  several  Ionization 
constants  are  for  horse  hemoglobin  at  25°C.  Jf  similar  data 
for  the  ionization  of  human  hemoglobin  becomes  available.  It 
can  be  Incorporated  with  possibly  some  Improvement  In  the 
representation  of  "acidity."  In  Model  II  (see  Appendix  D) , 
we  have  followed  an  approach  suggested  by  Pauling^ with 
regard  to  the  Ionization  of  the  hemoglobin  molecule. 

The  "elements"  In  rows  10  and  11,  protein  In  plasma  and 
In  red  cells  (other  than  hemoglobin),  are  Inserted  In  the 
model  at  this  stage  In  terms  of  moles  of  average  molecular 
weight  (equivalent  moles)  to  represent  the  wide  range  of 
molecular  species  of  protein  that  actually  occur  In  the  blood. 
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At  a  later  stage,  and  when  the  data  permit,  these  protein 
species  may  be  Introduced  separately  and  their  Ionization  may 
be  represented  In  the  same  fashion  as  suggested  for  the  case 
of  hemoglobin. 

As  mentioned,  row  12  contains  the  charge  equation,  which 
represents  the  conservation  of  charge  and  accounts  for  the 
Qlbbs-Donnan  equilibrium  across  the  red— cell  membrane.  At 
present  we  have  taken  b^2  0,  which  dictates  electrical 

neutrality  among  the  several  compartments.  Any  charge  dif¬ 
ference  that  may  occur,  as  between  red  cells  and  plasma,  can 
be  accounted  for  and  Its  Influence  on  relative  composition 
represented  by  the  Introduction  of  the  appropriate  molal 
charge-difference  value  of  b^^^*  Charge  differences  of  the 
order  of  — l6.8  millivolts  between  red  cells  and  plasma  are 
reported. Charge  differences  of  this  magnitude  could  have 
a  measurable  and  possibly  physiologically  significant  effect 
on  the  composition  of  plasma  and  red  cells.  Exploration  of 
this  phenomenon  by  means  of  the  model  remains  for  the  future. 
With  b^2  ■“  chloride  shift*  occurs,  so  that 

L°^Tpla37^  _  eella 

plasma  '  cells 


*It  can  be  shown  that  this  same  relation  holds  If  we 
allow  the  net  charge  b^2  Compartment  II  to  be  variable  and 

assiune  that  the  energy  required  to  maintain  this  charge  can 
be  expressed  by  adding  a  new  term  to  the  free-energy  function. 

1'.; 
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Instead  of  the  separate  relations 


plasma 

l-^^~3red  cells  * 

and 

[ncop  plasma 

^H^O^red  cells  * 

which  would  hold  if  the  charge 

equation  were  removed. 

Tlie  input  values  of  Vj^,  . 

*  *  ‘  ^12*  ^1#  •  •  •  > 

and  ,  bj^2»  which  the  equilibritun— solution  values 

,  x^Q  vfere  obtained,  are  listed  In  Table  2  In  the 


text 
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Appendix  C 

DETERMINATION  OF  FREE-ENERGY  VALUES 


So  far,  we  have  described  the  operation  of  the  model  In 
terras  of  the  several  compartments,  the  thermodynamics  Involved, 
the  Inputs,  and  the  mass  and  charge  conservation  constraints. 

We  shall  now  discuss  the  numerous  chemical  reactions  Involved, 
how  they  are  represented  In  the  model,  and  how  the  mass-action 
or  partlal-molal  free— energy  constraint  operates  (with  the 
others)  to  yield  an  equilibrium  solution.  The  mass— action 
equations  describe  the  formation  of  the  molecular  species  of 
columns  x^,  .  .  .  ,  x^q  In  the  three  compartments,  I:  Air  Out, 
II;  Arterial  Plasma,  and  III:  Red  Cells,  and  relate  similar 
species  occurring  In  more  than  one  compartment. 

For  example,  the  entry  1  In  column  12,  row  4,  Compartment 
II,  the  entry  1  In  the  same  column,  row  5,  and  the  entry  c^g 
In  the  bottom  row  together  express  the  chemical  equation, 

H"*’  +  OH"  ^  H2O  , 

In  the  mass— action  form, 

[OH"]  II 
-III 

where 


In 


[H+] 


II 


[He 


c 12  —  In  , 

and  Is  the  Ionization  constant  for  water  (on  the  mole- 
fraction  scale)  at  57°C.  It  was  shown  previously  that 
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constants  like  derived  from  equlllbrliun  constants  may  be 
used  In  the  model  with  the  seime  results  as  If  the  Individual 
P°/(RT)  values  for  each  species  were  obtained  from  tables  and 
used. 

In  a  manner  similar  to  that  described  for  the  water  re¬ 
actions  In  Appendix  A,  the  other  reactions  Involved  are 
entered  In  the  model.  For  the  same  species  occurring  In  both 
Compartments  II  and  III,  the  assumption  Is  made  that  their 
original  energy  coefficients  Cj  (not  necessarily  "Cj  values) 
are  the  same  In  both  compartments.  This  results,  for  example. 
In 

t^a^II  “  ^^°2^III* 

The  other  reactions  and  their  Cj  values  are  given  In  Table  4. 
In  this  table,  the  species  are  shown  In  square  brackets,  [  ], 
Indicating  concentration  In  mole  fraction.  The  Cj  values  were 
calculated  from  the  K  values  shown  In  Table  5,  all  converted 
to  the  mole— fraction  scale. 

It  Is  Important,  in  applying  these  data  to  the  model, 
that  the  K  values  from  the  literature  all  be  converted  to  the 
same  scale  (l.e.,  molar,  molal,  or  mole  fraction).  Most  K 
values  are  reported  on  the  molal  scale,  some  on  the  molar; 
and  In  some  cases,  unfortunately,  the  scale  Is  not  Indicated. 
All  K  values  were  converted  to  the  mole— fraction  scale  for 
use  In  the  model.  It  Is  also  Important  that  the  K  values  are 
those  for  the  temperature  of  Interest  (510.18  Kelvin  In  the 


Q?able  4 


RM-2519 

77 


REACTIONS  AND  FREE  ENERGY  VALUES^  FOR  MODEL  I,  57°C 


Column'^ 

I-Iass  Action 

Cj  a  In  K 

(1) 

In  [02]ii  / 
or  III 

-10.89 

(2) 

In 

or  III 

-7.69 

(5) 

In  [|N23  J--J-  /  [Ng]  j 
or  III 

-11.49 

(12)(24) 

.  In  j  j  COH~]  •j-.j-  /  TT 

or  III  or  III  or  III 

-59.23 

(4) 

In  [ohIjj  /  [HgOjj 

or  III  or  III 

=  c^2  +  in  [HgOljj  /  [HgOlj 
or  III 

-36.44 

(15)(25) 

In  [COgliT  [OHIt-j  /  [HC0-]  j 
or  III  or  III  or  III 

-21.20 

(29) 

in  [02]iii  [hBIIjii  /  [HBlOjjjj 

“  c^  +  In  [Og]  j  ^HBl  3  j  /  l-^®10^  j j j 

-16.23 

(i4)(26) 

In  [COglii  /  tH2C0^]jj 

or  III  or  III  or  III 

0 

(15) (27) 

In  [CO  ]  [OHl  /  [C0=]j  [h+] 

or  III  or  III  or  III  or  III 

6.25 

®A11  Cj  values 

to  In  Appendix  A. 


3hovm  are  the  a<ijusted 


values  referi*ed 


‘^The  values  of  Cj  =>  0  for  columns  (3)^  (6), 

(9).  (10),  (11),  (16),  (17),  (18),  (19),  (20), 
(23),  (28),  (30).  .  , 


(7),  (8), 
(21),  (22), 
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Table  5 


EQUILIBRIUM  CONSTANTS® 

Column 

Equilibrium  Constant 

(1) 

(2) 

(5) 

(12) (24) 

K  «  0.0232'’ 

K  -  0.5672*’ 

K  -  0.0127" 
logio  K  -  -13.55^ 

(^) 

logio  “  ^°®io  'v  '°®io  ^''2®^  II  / 

or  III 

logio  t''2°3ii  /  CH2°3i  "  1-21596® 
or  III 

(15) (25) 

log^O  ^  “  ^°8i0  ^  ~  '^‘’^lO  '^Ijj  pfj  * 

2 ' 

log-iQ  K.  ™  6.09 

•^H2C0j 

(29) 

logj^Q  K  -  -7-0486® 

(14)(26) 

log^Q  K  -  700" 

(15)(27) 

logio  "  '^^^lO  f^^°2pll 

H2CO3 

log-o  Kp  -  -10.25^ 

^COj 

®A11  K  values  on  the  molar  or  molal  scale  were  converted  to 
the  mole  fraction  scale  for  calculating  c*s. 

'’Ref.  3,  p.  3^- 

®Ref.  l6,  p,  1092. 

‘^Ref.  17,  p.  152. 

®Ref.  16,  p.  1465. 

^Ref.  9,  P-  272. 

®Ref.  3,  p.  64. 

"Ref.  3,  p.  80. 

^Ref.  16,  p.  1198. 
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model  30  far) ,  or  are  converted  to  the  value  for  .18 

the  necessary  thermodynamic  data  permit.  At  a  later  stage, 
the  appropriate  thermodynamic  functions  can  be  Introduced  Into 
the  model,  relating  free  energy  to  temperature  and  pressure, 
and  perhaps  to  electrical  and  gravitational  fields  and  to 
surface  effects. 


Appendix  D 
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MODEL  II  —  RESPIRATORY  MODEL  WITH  HEMOGLOBIN  STRUCTURE 

Mention  was  made  in  Sec.  VII  of  certain  deficiencies  in  our 
Model  I  of  the  respiratory  function.  Of  those  relating  to  the 
chemistry  of  the  blood,  the  most  Important  are  the  representa¬ 
tion  of  hemoglobin,  including  the  carbamino  reaction,  and  the 
multistaged  ionization  of  protein  and  protein— like  substances 
in  the  plasma  and  red  cells.  An  expanded  Model  II,  shown  in 
Table  6,  was  constructed  to  demonstrate  one  of  several  ways  in 

* 

which  the  repre;:entation  of  these  phenomena  might  be  improved. 

In  this  model,  two-staged  ionization  of  proteins  is  ac¬ 
counted  for  in  both  plasma  and  red  cells.  The  representation 
of  this  type  of  phenomenon  could  be  expanded  to  any  degree  of 
detail  desired.  In  addition  to  proteins,  it  might  Include  a 
wide  range  of  individual  molecular  species  such  as  amino  acids; 
and  it  might  Include  other  reactions  than  ionization  —  e.g., 
esterification. 

In  Model  II,  the  more  detailed  structuring  of  hemoglobin 
and  its  reactions  is  represented  by  a  separate  compartment  IV, 
as  if  a  separate  phase  were  involved.  There  are  distinct  ad¬ 
vantages  to  this  form  of  representation  for  such  complex 
structures  and . reactions  whenever  it  can  be  assumed  that  re¬ 
actions  with  different  parts  of  the  molecule  are  operating 

« 

A  number  of  hypotheses  exist  in  the  literature;  we  have 
chosen  those  due  to  Adair  (5)»  Pauling  (15)»  and  Roughton  (5)» 
among  others,  for  Illustrative  purposes  only.  In  cases  where 
there  are  competing  hypotheses,  alternative  models  could  be  set 
up  and  used  to  see  which  best  predicts  the  complex  obseirved 
phenomena. 


HM-2519 


RM-2519 

85 


independently.*  Formal  compartmentalizatlon  is  not 

essential  —  by  a  somev;hat  different  arrangement,  hemoglobin 

and  its  reactions  could  be  Incorporated  in  Compartment  III, 

Red  Cells.  Columns  5^—57  and  v;ere  designated  as  species 

of  hemoglobin,  and  as  portions  of  its  molecule,  to  correspond 

with  Adair's  hypothesls^^^  of  four  levels  of  oxidation  of  the 

hemoglobin  molecule  and  vilth  Linus  Pauling's  square— heme- 

structure  hypothesis  and  techniques  for  accounting  for  the 

variable  affinity  of  oxygen  with  changes  in  oxygen  concentra— 

( Ih) 

tion  and  v/ith  changes  in  pH.'  ' 

Coi-umns  58,  59,  and  54  account  for  the  CO^-^arbamino  link¬ 
ing  with  hemoglobin  and  its  relation  with  oxygen  and  hydrogen- 
ion  concentration.  The  mechanisms  suggested  by  Houghton 
(pp.  86—88  of  Ref.  4)  have  been  followed  in  this  representation. 

These  more  detailed  reactions  of  hemoglobin  as  repre¬ 
sented  in  compartment  IV  may  be  described  briefly  as  follovfs: 
The  six  possible  molecular  arrangements  of  hemoglobin  accord¬ 
ing  to  Pauling  are  here  designated  as  "B",  02"B",  0j^"B2"(a), 
0jj^"B2"  (t>) ,  Og"B^",  and  Og"Bij"  (columns  59— For  those 
molecular  arrangements  (columns  4l,  45,  44)  where  the  square 
oxygen— linked  heme  groups  Interact  nonsymetrically,  the  free 

*Wlthout  such  a  device,  it  vfould  be  necessary  to  set  up 
separate  columns  to  represent  perhaps  1  million  distinct 
variants  of  the  hemoglobin  molecule  corresponding  to  the 
different  ways  its  more  than  100  groups  and  parts  can  ionize, 
or  bond  with  COg  or  Og. 
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energy  l3  decreased  algebraically  by  a  factof,  RT  In  a,  for 
each  unbalanced  side,  v^here  a  Is  an  empirical  value  determined 

I 

from  experimental  evidence.  Inasmuch  as  AF  —  RT  In  K  , 

I 

where  K  Is  the  free— energy  change  accompanying  the  addition 
of  oxygen  to  heme,  the  "c"  values  for  the  model  may  be  ob¬ 
tained  from  the  logarithms  of  the  separate  and  appropriate  a 
factors  for  each  of  the  molecular  arrangements. 

In  an  analogous  manner,  Pauling  suggests  that  the  effects 
of  changes  In  hydrogen— Ion  concentration  on  hemoglobin  may  be 

I 

effected  by  an  empirical  correction  to  the  K  value.  In  this 
case,  empirical  factors  RT  In  3  and  a  so-called  acid— strength 
constant,  RT  In  A,  are  used  to  account  for  the  difference  In 
hydrogen— Ion  dissociation  of  oxygenated  and  unoxygenated 
hemoglobin.  These  correction  factors  are  determln'^d  from 
experimental  evidence.  The  precise,  mechanism  for  doing  this 
Is  discussed  below. 

In  terms  of  Model  II,  columns  54  and  55,  "A~"  and  H"A", 
represent  the  groups  of  the  hemoglobin  molecule  that  can 
Ionise  and  that  are  not  affected  by  oxygen  bonding.  There 
are  n^  of  these  hydrogen— bonding  groups.  The  energy  constant 
per  group  Is  assumed  to  be  Pauling's  In  A.  Coliimns  56  and 
57»  "C"””  and  H"C",  represent  an  additional  eight  Ionizing 
groups  of  the  hemoglobin  molecule  near  to  and  affected  by 
oxygen  bonding  ( "oxylablle") .  Ti\’o  acid  groups  Interacting 
with  each  heme  v;ere  postulated  by  Pauling  to  obtain  satis¬ 
factory  agreement  between  the  model  and  experimental  results. 
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Each  "oxylablle"  group  Is  assumed  also  to  have  the  energy 
constant  In  A  v/hen  no  oxygen  is  attached  on  the  heme  and 
In  A  -  In  p  when  oxygen  is  present.  Following  Roughton,  we 
allow  one  COg  to  attach  to  one  of  two  oxylabile  groups  If  no 
Og  is  present. 

I 

From  the  empirically  determined  values  for  K  ,  a.  A,  and 
3  from  Pauling,  it  is  therefore  possible  to  calculate  the 
values  c^jj^  —  C2^2|  Compartment  iv  of  Kodel  II.  The 

mass  balance  of  all  the  hemoglobin-associated  species  of 
Compartment  IV  is  related  to  total  hemoglobin  in  Compartment 
III,  Red  Cells,  through  column  35  of  Red  Cells  and  through 
the  Og,  COg,  and  H'*’  coefficient  entries  of  the  hemoglobin- 
associated  Input  "Elements." 

Model  II  has  not  been  validated;  until  this  is  done  it 

I 

cannot  be  certain  that  K  ,  a.  A,  and  p  values  have  been 
correctly  Interpreted  in  the  model.  For  example,  the  energy 
constants  Cj  of  Model  II  are  on  a  per-group  basis,  whereas 
Pauling's  A  and  3  values  may  apply  to  the  computation  of  a 


pair  of  groups 


RM-2519 


87 


Appendix  E 

SIMULATING  THE  RESPIRATORY  FUNCTION  ON  AN  ANALOG  COMPUTER 
THE  USES  OF  THE  COMPUIERS 

Computer  solutions  of  the  lung  model  have  been  undertaken 
with  two  goals  In  mind.  The  first  Is  the  obvious  necessity 
to  make  sensitive  checks  on  the  mathematical  model.  A 
physiological  model  Is  necessarily  founded  on  many  parameters 
and  a  number  of  hypotheses,  both  explicit  and  Implicit.  The 
ability  to  make  fast,  deliberate  tests  of  these  underlying 
assiimptlons  has  aided  materially  In  the  validation  of  this 
first  model. 

The  second  goal  Is  to  obtain.  In  physical  hardware,  a 
simulation  of  the  time-dependent  system.  Such  a  simulation 
should  present  the  complex  Interplay  of  component  chemical 
reactions;  but  It  should  do  more.  It  should  be  amenable  to 
the  analysis  of  Individual  reactions  as  v;ell  as  the  whole; 

It  should,  of  course,  compute  the  equillbrlian  conditions;  and 
It  should  have  sufficient  capacity  and  speed  to  enable  an 
analysis  of  the  transient  change  from  one  equlllbrlian  state 
to  another.  Such  a  change  In  equilibrium  conditions  will 
occur,  for  example,  when  the  lungs  are  suddenly  exposed  to  pure 
oxygen  Instead  of  the  normal  air  mixture. 

Almost  any  physiological  model  that  approaches  reality 
will  be  sufficiently  complex  to  obfuscate  one's  mental  Image 
of  the  system.  But  the  analog-computer  system  has  the 


ft 
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advantage  that  it  may  be  built  up  a  piece  at  a  time  to  any 
desired  complexity.  Hence,  the  goals  mentioned  above  seem 
perfectly  feasible,  although  certain  limitations,  such  as  com¬ 
puter  size,  are  apparent.  Despite  the  fact  that  the  two  goals 
overlap,  it  was  decided  to  attain  them  on  different  computers. 
The  digital  computer,  with  high  accuracy  but  relatively  long 
solution  time,  is  used  to  compute  the  equilibrium  conditions 
of  the  system  and  hence  to  test  the  underlying  hypothesis  of 
the  model.  The  analog  computer,  with  lower  accuracy  but  much 
shorter  solution  time,  can  be  used  to  search  through  a  range 
of  parameters  or  can  eventually  be  used,  we  believe,  to  simu¬ 
late  the  dynamic  aspects  of  the  model.  For  the  present,  it 
is  a  valuable  tool  for  tracing  the  effects  of  continuously 
varying  Inputs  to  the  system.  For  example,  the  analog  can  be 
used  in  complex  experimental  situations  for  estimating  equi¬ 
librium  constants  by  regarding  these  constants  as  parameters 
to  be  varied  continuously  on  the  machine  until  a  good  fit  is 
found  between  observations  and  predictions. 

CONTINUOUS  OPTIMIZATION  AND  THE  ANALOG  COMPUTER 

The  problem  of  continuous  optimization,  as  compared  to 
the  static,  one-time,  optimization  of  Appendix  D,  is  to  com¬ 
pute  and  display  the  transient  phenomena  of  a  complex  system 
ns  it  changes  under  the  Influence  of  a  time— dependent 
parameter.  The  central  postulate  is  Hamilton's  principle, 
the  principle  of  least  action  from  classical  mechanics. 
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The  system  Is  defined  as  In  the  text.  The  free— energy 
function  P(x)  Is  to  be  minimized  subject  to  mass— balance  and 


charge  constraints.  Using  the  Lagranglan,  we  write,  for  k  =  I, 
II.  Ill, 

ra  n  X.  m 

Q(x,ir)  =  F(X)  -  Z  TT.  <}).  ■  Z  X,  (c.+ln  :::i)  -  Z  ir.  (x),  (24) 

1=1  ^  1=1  J  J  X.  1=1  1  ^ 


n 


where  ({),(x)  =  —  b^.  At  equilibrium,  we  require 


-  c.  +  In  3:^  —  Z  a.  .T.  “  0  lor  each  J, 
dxj  J  Xj^  1=1  ^ 


(25) 


=  (J).  ( X)  ■=  0  for  each  1 , 


where  the  are  the  Lagranglan  multipliers. 

A  unique  solution,  or  equilibrium  condition,  obtains  on 
the  computer  when  a  set  of  mole  numbers  X  and  a  vector  r 
Identically  satisfy  the  equations  (25). 

If  the  physical  system  represented  by  Model  I  Is  well 
behaved.  It  will  exhibit  at  a  solution  certain  stability 
properties,  and  the  system  will  remain  stationary  until  a 
parameter  of  the  system  Is  redefined,  e.g.,  the  quantity  of 
a  reactant.  At  such  a  time,  transient  phenomena  are  exhibited 
until  a  new  equlllbrlvim  obtains.  For  most  parameters  of  the 
system,  small  changes  of  the  parameter  reflect  small  changes 
In  the  solution;  that  Is,  the  solution  Is  a  continuous  func¬ 
tion  of  the  parameters  of  variation.  Thus,  If  a  parameter  Is 
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a  slowly  changing  function  of  time,  the  solution  moves  In  a 
trajectory  that  continuously  minimizes  F(x).  Thus  (25)  always 
holds,  although  G  =  G(X,ir,t). 

THE  TECHNIQUE  OF  THE  ANALOG  SOLUTION 

Accurate  continuous  simulation  of  a  chemical  complex 
may  be  accomplished  on  an  analog  computer  by  means, 
essentially,  of  the  method  of  steepest  descent.  The  analog 
Is  desirable  because  of  the  characteristically  short  solution 
time  available,  particularly  In  the  high— gain  mode  of  compu¬ 
tation.  The  general  method  Is  shown  In  a  previous  pai>er, 
"Continuous  Programming  Methods  on  the  Analog, which 
contains  a  bibliography.  Some  modifications  have  ensued, 
however,  which  represent  a  departure  from  previous  pro¬ 
gramming  methods;  some  details  will  accordingly  be  shown  In 
this  Appendix. 

The  essential  Justification  for  the  method  employed  comes 
from  the  studies  of  convergence  questions  In  the  litera¬ 
ture,  ^  T.  Kose,  In  a  1956  paper,  considered  the 

convergence  of  a  solution  of  the  mathematical  programming 
problem  as  solved  on  an  analog  computer.  To  use 
his  result.  It  Is  necessary  to  postulate  the  existence  and 
uniqueness  of  a  solution,  considerations  which  are  taken  up 
for  more  general  programming  problems  elsewhere  In  the 
literature  of  mathematical  programming.  With  these 
hypotheses,  we  have  the  following  result. 
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Thoorom.  ^  0(X,Tr)  Is  a  strictly  ooncavo,  differentiable 
function  of  X  as  v/ell  as  a  strictly  convex  and  differentiable 


function  of  ir  for  Xj  y  0,  then  all  solutions  of  the  system  of 
of  differential  equations. 


9 


1  f  X .  =  0  and  <  0  , 


If  TT,  =0  and  ^  0  , 
—  1  -  ' 


0 


othcnclse  , 


1 


a  =  a  positive  time  constant. 


necessarily  converge  to  a  unique  saddle  point# 

It  can  be  shovm  that  the  conditions  of  the  theorem  are 
satisfied  by  the  function  G.  Tne  constant  a  Is  arbitrary  and 
we  wish  to  make  It  very  small  so  that  the  solution  may  be 
obtained  in  a  short  time.  Of  course,  various  practical  con¬ 
siderations  obtain  on  the  analog  computer,  because  the  com¬ 
puter  Is  a  physical  device  with  characteristics  (e.g., 
frequency  llmitatlors)  of  its  ovm,  and  these  characteristics 
supervene  on  the  theoretical  mathematical  statement.  V/lth 
these  considerations  in  mind,  v/e  may  proceed  to  program  the 
problem. 
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It  Is  required  to  maintain  the  relations  (25)  Identically, 
Independent  of  a  new  time  variable  t.  To- this  end,  the  rela¬ 
tions  {25)  are  programmed  on  amplifiers  In  the  hlgh-galn  mode, 
the  outputs  of  which  are  the  vectors  X  and  ir.  In  order  to 
explain  this  mode,  we  start  with  the  basic  operational  unit 
of  the  computer,  the  hlgh-galn  amplifier,  represented  by  the 
symbol 

®1 

— >  Oq  =  -  n(ej  +  eg)  (26) 

where  p.  Is  the  "gain"  of  the  amplifier  and  Is  usually  of  the 
order  of  10^.  Note  the  sign  change  and  the  addition.  Since 
the  output  variable,  e^.  Is  limited,  for  practical  electronic 
reasons,  to  +  100  units,  the  sum  of  the  Inputs  Is  bounded  by 

J^ 

+  10  units;  1-e.,  the  sum  of  the  Inputs  must  be  very  near 
zero,  and  we  regard  this  device  as  a  restriction  on  the  In¬ 
puts  . 

In  order  to  control  the  sum  of  the  Inputs,  the  amplifier 
Is  normally  used  with  a  "negative  feedback"  loop.  The  design 
of  this  loop  defines  the  mathematical  operation  of  the  unit. 
Thus,  Just  as  an  example.  It  Is  possible  to  take  the 
(positive)  square  root  of  an  Input,  e^^,  by  Inserting  a 
"multiplier"  In  the  feedback: 
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->  Sq  =  -X  (27) 


«  +x 


A  sign  change  occurs  In  the  multiplier.  A  heuristic  explana¬ 
tion  of  the  operation  may  be  had  by  regarding  the  amplifier 
as  an  Integrator.  The  output,  l.e.,  the  Integral,  will  be 
steady  v^hen  the  sum  of  the  Inputs,  l.e.,  the  derivative.  Is 
zero.  Thus,  a(— x  +x  )  Is  a  differential  operator  for  e^, 
the  whole  system  Is  a  function  of  time,  and  as  t  gets  large 
a(— X  +x  )  goes  to  zero.  This  system  Is  stable  since  If  e^^ 

Is  given  a  positive  Increment,  e^  takes  a  negative  Increment 
2 

and  —X  decreases  bringing  the  simi  of  the  Inputs  again  to 
zero.  The  whole  operation  should  be  regarded,  now,  as 
determining  e^.  If  the  gain  Is  sufficiently  high,  say  10^, 
the/accuracy  is  very  good  since  small  errors  In  e^  cause 
very  significant  errors  In  the  zero  sum  of  tlie  Inputs  and 
In  the  right  direction  to  correct  e^. 

An  Important  thing  to  note  Is  that  e^  was  arbitrarily 
defined;  it  has  no  necessary  relation  to  e^^.  But  once  It  Is 
defined,  e^  =  f(ej^).  It  Is  necessary  to  constiTJct  the  Inverse 
function  generator,  e^^  =  f~^(eQ),  for  the  feedback  loop. 

Also  the  feedback  loop  may  Itself  be  Implicitly  defined. 

Using  such  self— regulating  loops,  we  may  Instrument 
the  complex  equlllbrlimi  problem,  and  the  computer  will 


94 


automatically  seek  the  equilibrium  values,  l.e.,  values  for 
vfhlch  each  of  the  Interconnected  loops  Is  satisfied.  Such  a 
solution  requires  a  few  milliseconds.  Thus,  we  Instrument 
the  mass— balance  equations. 


n 

a.  .X .  —  b.  •  0  , 
J=1  J  ^ 


and  name  the  outputs  ir 


1* 


(2a) 


The  first  of  equations  (2)  are  Instrumented  and  the  outputs 
named  +Xj ; 


(29) 


Note  that  It  Is  necessary  to  generate  In  Xj. 

Since  the  outputs  of  (28)  and  (29)  are  In  the  Inputs  of  (29) 
and  (28),  respectively.  It  Is  sufficient  to  multiply  by  the 
appropriate  coefficients  and  Interconnect;  each  Is  the  feed¬ 
back  for  the  other.  The  parameters  b^  and  Cj  are  available 
as  constants  or  as  functions  of  time. 

As  an  example  of  the  above  discussion,  the  Ng  (nitrogen) 
distribution  lends  Itself  well  since  It  Is  Interconnected 
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with  the  other  species  only  through  Xj^,  which  we  may  consider 
constant  for  the  present. 

With  reference  to  Model  I,  the  following  diagram  was 
constructed: 


(30) 


The  generators  for  In  Xj  are  not  shown, 
c^^  satisfy 


c 


7 


The  constants  Cj  and 


for  the  reasons  given  In  the  text;  automatically  takes  the 
proper  sign;  +  ln(x^/  x^.)  is  a  constant  since  the  concentra¬ 
tion  of  Ng  in  alveolar  sac  air  is  a  postulated  constant  and 
c^  Is  the  tabulated  solubility  constant  for  Ng-  This  closed 
loop  Immediately  cof^putes  the  equilibrium  Xj  and  At 

equilibrium,  the  sum  of  the  Inputs  to  each  amplifier  Is  zero, 
as  required  of  amplifiers  In  the  hlgh-galn  mode,  and  hence 
equation  (25)  Is  satisfied.  Equations  Xj  ^  0  may  be  Invoked, 
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as  required,  by  means  of  a  limiter  around  the  Xj  amplifiers 
Actually  Xj  )>  0^  and  the  limiters  function  only  during  the 
turn— on  phase  of  the  computer.  Finally, 


X 


II 


^III 


or 


(51) 


In  — —  c,  ■»  In  K  , 

[X7]  ^ 


[X5] 


where  the  square  brackets  mean  "concentration"  and  the  mass- 
action  law  constant  is  satisfied. 


Computation  of  the  In  Xj  is  not  a  straightforward  pro¬ 


cedure.  It  is  done  here  by  a  linear  approximation,  but 
obviously  more  sophisticated  techniques  could  be  used.  Also, 
it  is  sometimes  necessary  to  attenuate  a  natural  high- 
frequency  phenomenon  that  occurs  in  high— gain  loops.  Very 
small  capacitors  may  be  used  without  affecting  the  computa¬ 
tion. 

The  basic  system  requires  Just  n  +  m  amplifiers;  but 

because  of  scale— factor  problems,  trace  elements  must  be 

computed  from  the  mass— action  constants  using  multipliers  as 

( 5) 

in  White,  et  al.,  or  by  some  other  means. 

THE  CONTINUOUS  ANALOG  SOLUTION 

Tlie  simplest  continuous  solution  for  the  equilibrium 


problem,  when  G  =  G(x,ir,  t),  can  be  Illustrated  by  means  of 
the  example  given  above. 
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Let 

=  b^(t)  =3  A  +  B  3ln  wt.  (^2) 

If  w  l3  a  moderate  frequency  —  say,  the  brcathtnc 
rate  —  then  the  computation  rate  of  the  machine  Is  so  high 
that  for  any  t  the  Instantaneous  equilibrium  conditions  will 
obtain.  The  function  Q  will  then  oscillate  with  angular 
frequency  w,  a  continuous  function  cf  b^(t).  Experiments  of 
this  sort  have  been  performed. 

Of  course,  each  of  the  parameters  Is  available  for  manual 
Investigation  of  Its  effect  on  the  system.  However,  the  more 
Interesting  cases  arise  when  one  attempts  to  define  the  delay 
functions,  diffusion  equations,  temperature  dependencies,  and 
such  complexities  as  may  actually  occur.  As  an  example, 
there  Is  the  question  of  the  equilibrium  trajectory  when  the 
pressure  Is  gradually  reduced  and  then  the  atmosphere  Is 
suddenly  switched  to  pure  oxygen.  Such  a  complex  question  has 
not  been  Investigated  but  seems  perfectly  feasible  of  solution. 

Methods  for  continuous  optimization  of  a  complex  system 
have  had  a  relatively  short  history  In  the  literature.  In 
general,  this  problem  Is  related  to  the  synthesis  of  optimal 
control  systems  or,  better,  the  optimal  synthesis  of  optimal 
control  systems.  In  such  a  problem,  a  complex  system  unde]>- 
golng  transformation  Is  automatically  to  be  optimized  with 
respect  to  a  quality  of  the  output.  Ordinarily,  then,  a 
feedback  link  Is  required  to  alter  one  or  more  Inputs,  and  a 
criterion  Is  required  to  measure  the  performance  of  the  system 
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v/ith  rcapect  to  a  (an  optimal)  standard.  In  particular,  see 
Kef.  21- 

The  present  problem  lo  somcv/hat  simpler  than  the  complete 
control  system,  bxit  Involves  more  complex  asjaocts  as  the 
physlolOGlcal  model  Increases  In  sl7,e.  For  the  present,  we 
continuously  (in  time)  compute  the  equilibrium  values  of  a 
complex  chemical  system  as  the  Inputs  vary  at  random.  later 
on,  the  Inputs  v/111  be  outputs  computed  from  other  physio¬ 
logical  subsystems,  and  the  whole  will  be  optimized  according 


to  viable  criteria. 


RKFERENCES 


nM-2519 

99 


1. 


Haldane,  J.  S. , 
Piaven,  Conn. , 


Respiration.  Yale  University  Press,  New 

1922. 


2. 


Best,  C.  H.,  and 
Holt  and  Co. , 


N.  B.  Taylor,  Tlie  Living  Body.  Henry 
Inc . ,  New  York^i  1944,  p  .  169. 


3.  Boothby,  Walter  M- ,  M.  D.  (ed.).  Handbook  of  Respiratory 
Physiology.  Air  University,  USAP  School  of  Aviation 
Medicine,  Randolph  Air  Force  Base,  September,  195^. 


4.  Roughton,  P.J.W.,  "Respiratory  Functions  of  the  Blood," 

In  Walter  H.  Boothby,  M.  D.  (ed.).  Handbook  of 
Respiratory  Physiology.  Air  Unlyerslty,  USAP  School 
of  Aviation  Medicine,  Randolph  Air  Force  Base, 

September,  195^ »  p.  76. 

5.  White,  W.  B. ,  S.  M.  Johnson,  and  Q.  B.  Dantzlg,  "Chemical 

Equilibrium  in  Complex  Mixtures,"  J.  Chem.  Phys . ,  Vol.  28, 
No.  5,  May,  1958,  pp.  751-755. 

6.  Brinkley,  S.  R. ,  "Note  on  the  Conditions  of  Equilibrium 

for  Systems  of  Many  Constituents,"  J.  Chem.  Phys., 

Vol.  14,  Nos.  9  and  11,  September  and  November,  19^6, 

PP •  563—564 ,  686 . 


7.  Krleger,  P.  J.,  and  W.  B.  VUilte,  "A  Simplified  Method 
for  Computing  the  Equlllbrliim  Composition  of  Gaseous 
Systems,"  J.  Chem.  Phvs . .  Vol.  I6,  No.  4,  April,  19^8, 
pp.  358-360. 


8.  Rossini,  Frederick  D.,  £t  ,  Selected  Values  of 

Chemical  Thermodynamic  Properties.  Circular  500, 
National  Bureau  of  Standards,  V/ashlngton  25^  D.  C., 
Febiniary  1,  1952. 

9.  Spector,  William  S.  (ed.).  Handbook  of  Biological  Data. 

Wright  Air  Development  Center,  USAP,  October,  19 56 . 


10.  Mac.  Innes,  D.  A.,  Donald  Belcher,  and  Theodore 

Shedlovsky,  "The  Meaning  and  Standardization  of  the 
pH  Scale,"  J.  Amer.  Chem.  Soc..  Vol.  60,  No.  5>  May, 

1938,  p.  1096^ 

11.  Latimer,  W.  M.,  and  J.  H.  Hildebrand,  Reference  Book 

of  Inorganic  Chemistry.  Macmillan  Company,  Inc., 

New  York,  1950,  p.  219. 


RM-2519 

100 


12.  Chilton,  A.  B.,  and  R.  B.  Stacy,  "A  Mathematical  Analysis 
of  Carbon  Dioxide  Respiration  In  Man,"  Bull.  Math. 
Blophvs. .  Vol.  14,  1952,  pp.  1-18. 

15.  Albritton,  Errett  C.  (ed.).  Standard  Values  In  Blood.  Air 

Force  Technical  Report  6059/  V/rlght  Air  Development 
Center,  USAF,  July,  I95I. 

l4.  Alberty,  R.  A.,  "The  Ionization  of  the  Hems— linked  Groups 
of  Hemoglobin,"  J.  Amer.  Chem.  Soc..  Vol.  77,  No.  17, 
September,  1955,' PP-  4522-4524. 

15-  Pauling,  Linus,  "The  Oxygen  Equilibrium  of  Hemoglobin  and 
Its  Structural  Interpretation,"  Froc.  Nat.  Acad.  Scl.. 
Vol.  21,  1955,  pp.  186-191. 

16.  Lange’s  Handbook  of  Chemistry.  9th  ed..  Handbook 

Publishers,  Inc.,  New  York,  1956,  p.  1465- 

17.  Washburn,  Edward  W.,  International  Critical  Tables. 

Vol.  6,  McGrav;— Hill  Book  Co.,  Inc.,  New  York,  1929, 
p.  152. 

18.  DeLand,  E.  C.,  Continuous  Programming  Methods  on  the 

Analog .  The  RAND  Corporation,  Paper  P— 1615, 

September  29,  1959* 

19.  Kose,  T. ,  "Solutions  of  Saddle  Value  Problems  bv 

Differential  Equations,"  Econometrlca .  Vol.  24,  1956, 

PP-  59-70. 

20.  Arrow,  K.  J. ,  and  L.  Hut^Icz,  A  Gradient  Method  for 

Approaching  Saddle  Points  01  Constrained  Maxima. 

Tlie  RAND  Corporation,  Paper  P-225,  June  I5,  1951- 

21.  Kolman,  R.  E.,  and  R.  W.  Koepske,  "Digital  Computers  In 

the  Optimization  of  Chemical  Reactors,"  Proc.  VJestern 
Joint  Computer  Conference.  Institute  of  Radio  Engineers, 
San  Francisco,  March,  1959,  P-  107- 


O'**' 


UNCLASSIFIED 


UNCLASSIFIED 


